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Abstract. We give an overview of the integrability of the Hirota-Kimura discretization 
method applied to algebraically completely integrable (a.c.i.) systems with quadratic 
vector fields. Along with the description of the basic mechanism of integrability (Hirota- 
Kimura bases) , we provide the reader with a fairly complete list of the currently available 
results for concrete a.c.i. systems. 
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1. Introduction 

The discretization method studied in this paper seems to be introduced in the geometric 
integration literature by W. Kahan in the unpublished notes [10J. It is applicable to any 
system of ordinary differential equations for x : K — » W 1 with a quadratic vector field: 

x = Q(x) + Bx + c, 

where each component of Q : lR n — > W 1 is a quadratic form, while B 6 Mat nxra and 
cGR" Kahan's discretization reads as 

- — - = Q(x,x) + )-B(x + x) + c, (1) 

where 

Q(x, x) = - (Q(x + x) - Q(x) - Q(x)) 

is the symmetric bilinear form corresponding to the quadratic form Q. Here and below 
we use the following notational convention which will allow us to omit a lot of indices: 
for a sequence x : Z — > R we write x for Xk and x for Xk+i- Eq. (TjQ) is linear with respect 
to x and therefore defines a rational map x = f(x,e). Clearly, this map approximates 
the time-e-shift along the solutions of the original differential system, so that Xk ~ x(ke). 
(Sometimes it will be more convenient to use 2e for the time step, in order to avoid 
appearance of various powers of 2 in numerous formulas.) Since eq. ([1]) remains invariant 
under the interchange x -H- x with the simultaneous sign inversion e i— > — e, one has the 
reversibility property 

f~ 1 {x ) e) = f(x, -e). 

In particular, the map / is birational. 

W. Kahan applied this discretization scheme to the famous Lotka-Volterra system and 
showed that in this case it possesses a very remarkable non-spiralling property. This 
example is briefly discussed in [17J. Some further applications of this discretization have 
been explored in 

The next, even more intriguing appearance of this discretization was in the two pa- 
pers by R. Hirota and K. Kimura who (being apparently unaware of the work by Kahan) 
applied it to two famous integrable system of classical mechanics, the Euler top and the 
Lagrange top PIT2]. For the purposes of the present text, integrability of a dynamical 
system is synonymous with the existence of a sufficient number of functionally indepen- 
dent conserved quantities, or integrals of motion, that is, functions constant along the 
orbits. We leave aside other aspects of the multi-faceted notion of integrability, such as 
Hamiltonian ones or explicit solutions. Surprisingly, the Kahan-Hirota-Kimura discretiza- 
tion scheme produced in both the Euler and the Lagrange cases of the rigid body motion 
integrable maps. Even more surprisingly, the mechanism which assures integrability in 
these two cases seems to be rather different from the majority of examples known in the 
area of integrable discretizations, and, more generally, integrable maps, cf. [25]. The case 
of the discrete time Euler top is relatively simple, and the proof of its integrability given 
in [S] is rather straightforward and easy to verify by hands. As it often happens, no ex- 
planation was given in [S] about how this result has been discovered. The "derivation" of 
integrals of motion for the discrete time Lagrange top in [12] is rather cryptic and almost 
uncomprehensible . 
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We use the term "Hirota-Kimura (HK) type discretization" for Kalian discretization in 
the context of integrable systems. At the Oberwolfach meeting on Geometric Integration 
in 2006, T. Ratiu proposed to apply the Hirota-Kimura discretization to the Clebsch case 
of the rigid body motion in an ideal fluid and to the Kovalevsky top. The claim on the 
integrability of the HK discretization of the Clebsch case was proven only several years 
later in [T7], while the integrability of the HK discretization of the Kovalevsky top remains 
an open problem (although there are some indications in favor of its non-integrability) . 
Anyway, the general question of integrability of the HK type discretizations turns out to 
be very intriguing and rather difficult. In the present overview, we will present a rather 
long list of examples of integrable HK discretizations. Actually, this list is so impressive 
that in [T7j we conjectured that HK discretizations of algebraically integrable systems 
always remain integrable. At present, we have some indications that this conjecture 
is wrong (although a rigorous proof of non-integrability for any of the apparently non- 
integrable cases remains elusive). Nevertheless, the sheer length of our list of examples 
clearly shows that there exist some general mechanisms that ensure integrability at least 
under certain additional assumptions. We think that to uncover general structures behind 
the integrability of HK discretizations is a big and important challenge for the modern 
theory of (algebraically) integrable systems. 

The structure of our overview is as follows. In Section [2] we demonstrate some general 
sufficient conditions for a HK discretization to be measure preserving. These conditions 
are not related to integrability and do not cover all special cases considered in the main 
text. All our examples turn out to be measure preserving but for the majority of them we 
only can prove this property individually and do not know any general mechanisms. In 
Section [3] we present a formalization of the HK mechanism from [T2], which will hopefully 
unveil its main idea and contribute towards demystifying at least some of its aspects. We 
introduce a notion of a "Hirota-Kimura basis" (HK basis) for a given map /. Such a basis 
$ is a set of simple (often monomial) functions, $ = (ipi, . . . , tpi), such that for every orbit 
{f l {x)} of the map / there is a certain linear combination C\(p\ + . . . + c\tp\ of functions 
from $ vanishing on this orbit. As explained in Sect. |3j this is a new mathematical 
notion, not reducible to that of integrals of motion, although closely related to the latter. 
We lay a theoretical fundament for the search for HK bases for a given discrete time 
system, and discuss some practical recipes and tricks for doing this. Sections HHTB1 contain 
our list of examples of algebraically integrable systems with quadratic vector fields, for 
which the HK discretization preserves integrability. For all the examples we provide the 
reader with a rather complete set of currently available results. The proofs are omitted 
almost everywhere. One of the reasons is that our investigations are based mainly on 
computer experiments, which are used both for discovery of new results and for their 
rigorous proof. For a given system, a search for HK bases can be done with the help of 
numerical experiments. If the search has been successful and a certain set of functions 
$ has been identified as a HK basis for a given map /, then numerical experiments can 
provide a very convincing evidence in favor of such a statement. A rigorous proof of such 
a statement turns out to be much more demanding. At present, we are not in possession 
of any theoretical proof strategies and are forced to verify the corresponding statements 
by means of symbolic computations which in some cases turn out to be hardly feasible due 
to complexity issues. All these issues are intentionally avoided in our presentation here; 
an interested reader may consult [T7] for a detailed discussion of one concrete example 
(Clebsch system). Our main goal here is to document the available results on integrable 
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HK discretizations and to attract attention of specialists in integrable systems and in 
algebraic geometry to these beautiful and mysterious objects which are definitely worth 
further investigation. 

2. Some general results on invariant measures 

In this section, we prove the existence of invariant measures for Kahan discretizations of 
two classes of dynamical systems with quadratic vector fields (not necessarily integrable). 
The first class reads: 

N 

x% = ciijxj + d, 1 < % < N, 

3=1 

with a skew-symmetric matrix A = (ay)f«- =1 = — A T . The Kalian's discretization reads: 

N 

Xi — Xi = e aijXjXj + eat, 1 < i < N. (2) 

3=1 

Proposition 1. The map x = f(x,e) defined by eqs. |IP has an invariant volume form: 

dx <f)(x, e) dx\ A ... A dx^ 

det ^~ = T( V ^ J u = w ' u = 77 \ ' 

ox (p{x,e) <P{x,£) 

where <p(x, e) = det(l — eAX), with X = diag(xi, . . . , xn), is an even polynomial in e. 
Proof. Equations 02]) can be put as 

x = A~ 1 (x,e)(x + ec), A(x, e) = 1 — eAX. 
As for any Kahan type discretization, there holds the formula 

dct dx = detA(x,-e) 
dx det A(x, e) 

Indeed, differentiation of equations (J2]) with respect to Xj will give the j-th column of the 
matrix equation 

Ox 

A ^,e)— = A(x, -e). 

It remains to notice that detA(x,e) = det A(x, — e). Indeed, due to the skew-symmetry 
of A, we have: det(l - eAX) = det(l - eX T A T ) = det(l + eXA) = det(l + eAX). □ 

The second class consists of equations of the Lotka-Volterra type: 

( r \ 



Xi X{ 



with a skew-symmetric matrix A = (ay)^\- =1 = — A T . The Kahan's discretization (with 
the stepsize 2e) reads: 

N 

Xi — Xi = ebi[x,i + Xi) + e 2J Oij{xiXj + XiXj), 1 < i < N. (4) 

3=1 
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Proposition 2. The map x = f(x,e) defined by equations 0) has an invariant volume 
form: 

dx X\X2 ■ ■ ■ xn ,* dxi A ... A 6xn 

det — = / lo = u, lo — . 

OX X\X2'--Xn X\X2'--Xn 

Proof. Equations @ are equivalent to 

00 q OCq 



1 + ebi + e Y,f=i a ij x j 1 — ebi — e J2jLi a ^ x • 



Vi- (5) 



We denote d^x, e) = 1 — efej — e ^,- =1 a^Xj. In the matrix form equation fll]) can be put 

as 



x = A _1 (x,e)(l + eB)x, 

where the i-th diagonal entry of A(x,e) equals di(x,e), while the (i,j)-th off-diagonal 
entry equals —ex^a^. In other words, A(x,e) = D(t — eYA), where D = D(x,e) = 
diag(di, . . . , djv) and Y = diag(yi, . . . y^). Formula ([3]) holds true also in the present 
case, and it implies: 

, dx det D(x, -e) det(l + eYA) 

QGt — 

dx detD(x,e) det(l-eFA)' 
The second factor equals 1 due to the skew-symmetry of A, while the first factor equals 

di(x, —e) ■ ■ ■ d N (x, —e) Xi---Xjq 
di(x, e) • • • dpf(x, e) Xi---x^ 

by virtue of (j^J). □ 

The statement of Theorem [2] for the Kahan discretization of the classical Lotka-Volterra 
model with two species, x = x(a — by), y = y(cx — d), was found in [23] and used to 
explain the non-spiralling behavior of the numerical orbits in this case. 

3. HlROTA-KlMURA BASES AND INTEGRALS 

In this section a general formulation of a remarkable mechanism will be given, which 
seems to be responsible for the integrability of the Hirota-Kimura type (or Kahan type) 
discretizations of algebraically completely integrable systems. This mechanism is so far 
not well understood, in fact at the moment we do not know what mathematical structures 
make it actually work. 

Throughout this section / : M. n — > M. n is a birational map, while hi, tpi : lR n — > K stand 
for rational, usually polynomial functions on the phase space. We start with recalling a 
well known definition. 

Definition 3. A function h : M n — > M. is called an integral, or a conserved quantity, 

of the map f , if for every iGl™ there holds 

h(f(x)) = h(x), 
so that h(f l (x)) = h(x) for all % £ Z. 

Thus, each orbit of the map / lies on a certain level set of its integral h. As a conse- 
quence, if one knows d functionally independent integrals hi, . . . , h^ of /, one can claim 
that each orbit of / is confined to an (n — <i)-dimensional invariant set, which is a common 
level set of the functions h\, . . . ,hd- 
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Definition 4. A set of functions $ = (tp\, . . . , (pi), linearly independent over M., is called 
a Hirota-Kimura basis (HK basis), if for every x G W 1 there exists a vector c = 
(ci, . . . , q) 7^ such that 

c m (f\x)) + ... + c m (f\x)) = (6) 

holds true for all i G Z. For a given x G M™, the vector space consisting of all c £ R 1 
with this property will be denoted by K$(x) and called the null-space of the basis $ (at the 
point x). 

Thus, for a HK basis $ and for c G K$(x) the function h = c\<p\ + ... + ci<pi vanishes 

along the /-orbit of x. Let us stress that we cannot claim that h = C\(p\ + + Ci(pi is 

an integral of motion, since vectors c G K$(x) do not have to belong to K<^(y) for initial 
points y not lying on the orbit of x. However, for any x the orbit {f l (x)} is confined to 
the common zero level set of d functions 

hj = c?V + . . . + cpVj = °> j = 1, . . . , <i, 

where the vectors = (c^ , . . . , c| ) G form a basis of K${x). We will say that the 
HK basis $ is regular, if the differentials dh\, . . . , dhd are lineraly independent along the 
the common zero level set of the functions hi, . . . , h$. Thus, knowledge of a regular HK 
basis with a d- dimensional null-space leads to a similar conclusion as knowledge of d inde- 
pendent integrals of /, namely to the conclusion that the orbits lie on (n — d)- dimensional 
invariant sets. Note, however, that a HK basis gives no immediate information on how 
these invariant sets foliate the phase space M n , since the vectors c^\ and therefore the 
functions hj, change from one initial point x to another. 

Although the notions of integrals and of HK bases cannot be immediately translated 
into one another, they turn out to be closely related. 

The simplest situation for a HK basis corresponds to I = 2, dim K^>(x) — d — 1. In 
this case we immediately see that h = <fi/<f2 is an integral of motion of the map /. 
Conversely, for any rational integral of motion h — (fi/(f2 its numerator and denominator 
<fi, <f2 satisfy 

c m (f(x)) + c 2lf2 (r(x)) = 0, iGZ, 

with c\ = 1, c 2 = —h(x), and thus build a HK basis with / = 2. Thus, the notion of a HK 
basis generalizes (for I > 3) the notion of integrals of motion. 

On the other hand, knowing a HK basis $ with dimK^,(x) = d > 1 allows one to find 
integrals of motion for the map /. Indeed, from Definition H] there follows immediately: 

Proposition 5. If $ is a HK basis for a map f , then 

K*(f(x)) =K*(x). 

Thus, the <i-dimensional null-space K$(x) G Gr(d,l), regarded as a function of the 
initial point x G M n , is constant along trajectories of the map /, i.e., it is a Gr(d,l)- 
valued integral. Its Pliicker coordinates are then scalar integrals: 

Corollary 6. Let $ be a HK basis for f with dim K<s>(x) = d for all i6R". Take a basis 
of K$(x) consisting of d vectors cW G M 1 and put them into the columns of a I x d matrix 
C(x). For any d-index a = (a>i, . . . , aj) C {1,2, ... ,n} let C a = C ai _ ad denote the d x d 
minor of the matrix C built from the rows ct\, . . . ,ctd- Then for any two d-indices a, (3 
the function C a /Cp is an integral of f . 



ON INTEGRABILITY OF HIROTA-KIMURA TYPE DISCRETIZATIONS 



7 



Especially simple is the situation when the null-space of a HK basis has dimension 

Corollary 7. Let $ be a HK basis for f with dimK$(x) = 1 for all x E W 1 . Let 

K$(x) = [ci(x) : . . . : q(x)] E MP' -1 . Then the functions Cj/ck are integrals of motion for 
/■ 

An interesting (and difficult) question is about the number of functionally independent 
integrals obtained from a given HK basis according to Corollaries E] and It is possible 
for a HK basis with a one-dimensional null-space to produce more than one independent 
integral. The first examples of this mechanism (with d = 1) were found in p2] and 
(somewhat implicitly) in [8]. 

We note, however, that HK bases appeared in a disguised form in the continuous time 
theory long ago. We mention here two relevant examples. 

• Classically, integration of a given system of ODEs in terms of elliptic functions 
started with the derivation of an equation of the type y 2 = P^y), where y is one 
of the components of the solution, and Pi{y) is a polynomial of degree 4 with 
constant coefficients (depending on parameters of the system and on its integrals 
of motion), see examples in Sections [121 HSJ This can be interpreted as the claim 
about $ = (y 2 , y 4 , y 3 , y 2 , y, 1) being a HK basis with a one-dimensional null- 
space. 

• According to [U Sect. 7.6.6], for any algebraically integrable system, one can choose 
projective coordinates y , yx, . . . , y n so that quadratic Wronskian equations are 
satisfied: 

n 

Ekl 

k,l=0 

with coefficients af!j depending on integrals of motion of the original system. 
Again, this admits an immediate interpretation in terms of HK bases consist- 
ing of the Wronskians and the quadratic monomials of the coordinate functions: 

®n = (inyj - ViVv {ykVi}k,i=o) ■ 

Thus, these HK bases consist not only of simple monomials, but include also more com- 
plicated functions composed of the vector field of the system at hand. We will encounter 
discrete counterparts of these HK bases, as well. 

At present, we cannot give any theoretical sufficient conditions for existence of a HK 
basis $ for a given map /, and the only way to find such a basis remains the experimental 
one. Definition 0] requires to verify condition ([6]) for all % E Z, which is, of course, 
impractical. However, it is enough to check this condition for a finite number of iterates 

Typically (that is, for general maps / and general monomial sets <&), the dimension of 
the vector space of solutions of the homogeneous system of s linear equations (jSJ) with 
i — io, io + 1, . . . , io + s — 1 decays with the growing s, from I — 1 for s — 1 down to 
for s — I. If, however, $ is a HK basis for / with a <i-dimensional null-space, then this 
dimension fails to drop starting with s — I — d + 1. Thus, the dimension of the solution 
space of the system (JBJ) with i = i Q , i + 1, . . . , i + s — 1 is equal to I — s for all 1 < s < I — d, 
and remains equal to d for s — I — d + 1. It is easy to see that this situation can be also 
characterized as follows: the d- dimensional spaces of solutions of the system ([6]) with 
i = io, io + 1, . . . , io + I — d — 1 and of the system ([6} with i = i Q + 1, i + 2, . . . , i + I — d 
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coincide with each other and with K$>(x). The most important particular case of this 
characterization corresponds to d — 1 and reads as follows: 

Proposition 8. A set $ = (tpx, . . . , ipi) is a HK basis for a map f with dimK<p(x) = 1 
if and only if the unique solution [c\ : . . . : q] G MP' -1 of the system with i = 
io, if) + 1, . . . , iq + I — 2 coincides with the unique solution of the analogous system with 
i — io + 1, io + 2, . . . , io + I — 1, in other words, if this unique solution [c\ : . . . : q] is an 
integral of motion: 

[c l (x):...:c l (x)] = [c 1 (f(x)):...:c l (f(x))]. 

At this point it should be mentioned that a numerical testing of the above criterium 
usually represents no problems, but the corresponding symbolic computations might be 
extremely complicated, due to the complexity of the iterates f % (x). While the expression 
for f(x) is typically of a moderate size, already the second iterate f 2 (x) becomes typically 
prohibitively big. See the detailed discussion of the complexity issue for the case of the 
Clebsch system in [I?]. In such a situation it becomes crucial to reduce the number of 
iterates involved in (Q as far as possible. Several tricks which can be used for this aim 
are also described in [17]. Here is one of them. 

Proposition 9. Consider the non-homogeneous system of I — 1 equations 

C lVl (f(x)) + ...+ Q_!^_!(f (x)) = Vl(f(x)), i = i Q ,i + 1, . . . ,i + I ~ 2. (7) 

Suppose that the index range % G [io, io + l — 2] in eq. |7|) contains but is non-symmetric. 
If the solution of this system (ci(x, e), . . . , q_i(x, e)) is unique and is even with respect to 
e, then all Ck(x, e) are conserved quantities of the map f , and $ = (<pi, . . . , cp{) is a HK 
basis for f with dimK$(x) = 1. 

Proof. Considering the non-homogeneous system ([7]) instead of the homogeneous one ([HD 
corresponds just to fixing an affine representative of the projective solution [c\ : . . . : q] 
by q = — 1. The reversibility of the map e) = f(x, — e) yields that equations of the 

system ([7]) are satisfied not only for i G [io, io + I — 2] but for i G [— (io + / — 2), —io), as 
well. Since, by condition, the intervals [io, io + I — 2] and [—(io + I — 2), — io] overlap but 
do not coincide, their union is an interval containing > / integers. □ 

Of course, it would be highly desirable to find some structures, like Lax representa- 
tion, bi-Hamiltonian structure, etc., which would allow one to check the conservation of 
integrals in a more clever way, but up to now no such structures have been found for any 
of the HK type discretizations. 

4. Weierstrass differential equation 

Consider the second-order differential equation 

x = 6x 2 — a. (8) 

Its general solution is given by the Weierstrass elliptic function p(i) = p(t,g2,g?) with 
the invariants gi = 2a, g^ arbitrary, and by its time shifts. Actually, the parameter g 3 
can be interpreted as the value of an integral of motion (conserved quantity) of system 

x 2 — 4x 3 + 2ax = —g-s- 
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Being re-written as a system of first-order equations with a quadratic vector field 

x = y, 
y = 6x 2 — a 



(9) 



equation (jHJ) becomes suitable for an application of the Kahan-Hirota-Kimura discretiza- 
tion: 

x-x= e -(y + y), (iq) 

y — y = e (6xx — a) . 

Eqs. ([ID]) , put as a linear system for (x,y), read: 

1 -e/2\ fx\ = fx + ey/2\ 
-6ex 1 J \y) \y-ea J' 

This can be immediately solved, thus yielding an explicit birational map (x, y) — f(x, y, e): 

x + ey — e 2 a/2 



x 



V 



1 — 3e 2 x 
y + e(6x 2 — a) + 3e 2 xy 



1 — 3e 2 x 

This map turns out to be integrable: it possesses an invariant two-form 

dxAdy 

u = T^V (12) 

and an integral of motion (conserved quantity): 

. y 2 — Ax 3 + 2ax + e 2 x(y 2 — 2ax) — e A a 2 x „. 

1 (x, V, e) = : ir- 2 • 13 

1 — 3e z x 

Both these objects are 0(e 2 )-perturbations of the corresponding objects for the continuous 
time system ([9]). The statement about the invariant 2- form (fl2"l) is not difficult to prove. 
The following argument exemplifies considerations which hold for an arbitrary Kahan 
discretization ([!]). Differentiating eqs. (JTUJ with respect to x and to y, we obtain the 
columns of the matrix equation 

1 -e/2\ d(x,y) _ ( 1 e/2 N 



whence _ 

d(x, y) 1 — 3e 2 x 



6ex 1 J d(x,y) \ 6ex 1 
det 



d(x, y) 1 — 3e 2 x 

This is equivalent to the preservation of ( I12p . The statement about the conserved quantity 
is most simply verified with any computer system for symbolic manipulations. 

System fflTJ]) is known in the literature on integrable maps, although in a somewhat 
different form. Indeed, it is equivalent to the second order difference equation 

x — 2x + x = e 2 (3x(x + x) — a) ^ x — 2x + x = — 



e 2 (6x 2 — <\ ) 



3e 2 x 

This equation belongs to class of integrable QRT systems [2T>P^]; in order to see this, one 



should re-write it as 

e 2 (6x 2 -a)(l + e 2 x) 

x - 2x + x =- 



1 - 2t 2 x - 3t A x 2 



10 MATTEO PETRERA, ANDREAS PFADLER, AND YURI B. SURIS 

This difference equation generates a map (x, x) y (x, x) which is symplectic, that is, 
preserves the two-form uj = dx A dx, and possesses a biquadratic integral of motion 

I(x, x, e) — (x — x) 2 — 2e 2 xx(x + x) + e 2 a(x + x) — e 4 (3x 2 x 2 — axx). 

Under the change of variables (x,x) h-> (x,y) given by the first equation in (iTTjl . these 
integrability attributes turn into the two-form (T12"]) and the conserved quantity (1131) (up 
to an additive constant). 

We note that a more usual QRT discretization of the Weierstrass second order equation 
(JED would be 

_ e 2 (6x 2 — a) , 

with a simpler conserved quantity 

J(x, x, e) = (x — x) 2 — 2e 2 xx(x + x) + e 2 a(x + x). 
Eq. ( 1T^|) is equivalent to 

x — 2x + x = e 2 (2x(x + x) + 2x 2 — a) , 
which is not obtained by the Kahan-Hirota-Kimura method. 

5. Some two-dimensional integrable systems 

5.1. The three-dimensional Suslov system. The three-dimensional nonholonomic 
Suslov problem [26] is defined by the following system of differential equations 

rh = m x u + Xa, (a,u) = 0, (15) 

where u = 002, W3) T is the angular velocity, m = (mi,m2,m 3 ) T = Iu is the angular 
momentum, / is the inertia operator, a is a unit vector fixed in the body, and A is the 
Lagrange multiplier. In a basis where a = (0, 0, 1) T and 



h 





h 3 





h 


I23 


/l3 


I23 


h 



the constraint (a, a;) = reduces to u 3 = 0, and equations of motion f TL3|) read 

hUi = -(/l3^1 + 

I 2 U 2 = (i"l3^1 + /23W2M, (16) 
+ = (h - h)^i^2 + A. 

The first two equations in (TIB]) form a closed system for oj\ and cj 2 - It possesses a 
conserved quantity H = Iiuf + After the solution of this system is found (Suslov 

gave it in terms of trigonometric and exponential functions), one finds the Lagrange 
multiplier A from the third equation in ffl6|) . 

To put the first two equations in f JT6|) into a more convenient form, one can introduce 
the coordinates x = I13UJ1 + /23W2, y = ^h^i — I13I2U2, and arrives at 

x = axy, , . 

2 17 

where a = \jl\l2- This system admits a conserved quantity 

H = x 2 + ay 2 . 
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Proposition 10. [5] HK discretization of system [Tfy , 

x — x = ea(xy + xy), 



y — y = —2exx, 

possesses an invariant two-form and an integral of motion, given by 

dx A dy x 2 + ay 2 

to = — — ; rr, Hie) 



x[x 2 + ay 2 ) 1 + e 2 ax 2 

Actually, the invariant two-form was not given in [5]; rather, this paper contains an 
explicit solution of the discrete time Suslov system and its qualitative analysis. 

5.2. Reduced Nahm equations. In [9] Nahm equations associated with symmetric 
monopoles are considered. Assuming rotational symmetry groups of regular polytops 
leads to solutions of Nahm equations in terms of elliptic functions. Reduced equations 
corresponding to tetrahedrally symmetric monopoles of charge 3, to octahedrally sym- 
metric monopoles of charge 4, and to icosahedrally symmetric monopoles of charge 6 
are two-dimensional algebraically integrable systems with quadratic vector fields. More 
concretely, the reductions of Nahm equations derived in [9] read: 

(i) tetrahedral symmetry: 

2 2 

x = x — y , . , 

y ' 18 
y = -2xy, 

with an integral of motion H = y(3x 2 — y 2 ); 

(ii) octahedral symmetry: 

x = 2x 2 - 12y 2 , . . 

y = -6xy - Ay , 

with an integral of motion H = y{2x + 3y)(x — y) 2 ; 

(iii) icosahedral symmetry: 

x = 2x 2 — y 2 , , . 

9 (20) 
y = -lOxy + y, 

with an integral of motion H = y(3x — y) 2 (Ax + y) 3 . 

HK discretizations of systems fll8H20p turn out to be algebraically integrable. 

Proposition 11. 

(i) HK discretization of system (TJSj], 

x — x = e(xx — yy), 



y-y= -e(xy + xy), 

possesses an invariant two-form and an integral of motion, given by 

dx Ady y(3x 2 - y 2 ) 

u = 77^73 H ( e ) 



y(3x 2 — y 2 ) ' 1 1 — e 2 (x 2 + y 2 ) ' 
(ii) HK discretization of system [W\) , 

x — x = e (2xx — V2yy) , 
y-y = -e(3xy + 3xy + Ayy) , 



12 



MATTEO PETRERA, ANDREAS PFADLER, AND YURI B. SURIS 



possesses an invariant two-form and an integral of motion given by 

dx A dy 



Hit) 



y(2x + 3y)(x - y)' 
y(2x + 3y)(x - y) 2 



1 - 10e 2 (x 2 + Ay 2 ) + e 4 (9x 4 + 272x 3 y - 352a;?/ 3 + 696y 4 ) 
(iii) HK discretization of system (G5J), 

x — x = e {2xx — yy) , 

y-y = t(-5xy-5xy + yy), 

possesses an invariant two-form and an integral of motion given by 

„ _ dx a dy Hie) - y( 3x - y) 2 ( Ax + y) 3 




y{3x — y)(Ax + y) ' 1 + e 2 C2 + e 4 C4 + e 6 C6 ' 

with 

c 2 = -7(5x 2 - y 2 ), 

c 4 = 7(37x 4 + 22x 2 y 2 - 2xy 3 + 2y 4 ), 

c 6 = -225x 6 + 3840x 5 y + 80xy 5 - 5Ux 3 y 3 - 19x 4 y 2 - 206x 2 y 4 . 

6. EULER TOP 

The differential equations of motion of the Euler top read 

CX1X2X3, 

a 2 x 3 xi, (21) 
a> 3 x 1 x 2 , 

with real parameters «j. This is one of the most famous integrable systems of the classical 
mechanics, with a big literature devoted to it. It can be explicitly integrated in terms of 
elliptic functions, and admits two functionally independent integrals of motion. Actually, 
a quadratic function H(x) = r )\x\ + 72X 2 + 73X3 is an integral for eqs. (j2Ti) as soon as 
7i«i + 72«2+72«2 = 0. In particular, the following three functions are integrals of motion: 



H\ = a 2 x 3 — «3X 2 , H 2 = «3X 1 — aix 3 , H 3 = ot\x 2 — ot2X 1 . 

Clearly, only two of them are functionally independent because of ot\Hi -\-a2H2 -\-a3H3 = 0. 
These integrals appear also on the right-hand sides of the quadratic (in this case even 
linear) expressions for the Wronskians of the coordinates xf 

x 2 x 3 - x 2 x 3 = Hxxi, 

x 3 xi - X3X1 = H 2 x 2 , (22) 
xix 2 - x\x 2 = H 3 x 3 . 

Moreover, one easily sees that the coordinates Xj satisfy the following differential equations 
with the coefficients depending on the integrals of motion: 

x\ = (H 3 + a 2 x 2 1 )(a 3 x 2 l - H 2 ), 
x\ = (Hi + a 3 a^)(aia^ - H 3 ), 
x\ = (H 2 + aixl)(a 2 xl - H x ). 
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The fact that the polynomials on the right-hand sides of these equations are of degree 
four implies that the solutions are given by elliptic functions. 
The HK discretization of the Euler top [8] is: 

= eai(x 2 x 3 + x 2 x 3 ), 

= ea 2 (x 3 x 1 + x 3 Xi), (23) 
= ea 3 (x 1 x 2 + x 1 x 2 ). 

(In this form it corresponds to the stepsize 2e rather than e.) The map / : x t— > x 
obtained by solving (T2151) for x is given by: 




x 




— eoc\X 3 —€OL\X 2 

f{x,e) = A~ 1 (x,e)x, A(x,e) = I — ea 2 x 3 1 — ea 2 Xi I . (24) 



2 \ 



-ea 3 Xi 1 I 

It might be instructive to have a look at the explicit formulas for this map: 

X\ + 2ea\X 2 x 3 + e 2 x\{— a 2 a 3 x\ + a 3 oi\x\ + ctia 2 x\) 



A(x,e) 




x 2 + 2ea 2 x 3 X\ + e 2 x 2 (a 2 a 3 x\ — a 3 oi\x\ 


+ aia 2 xl) 


A(x,e) 




x 3 + 2ea 3 X\X 2 + e 2 x 3 (a 2 a 3 x 2 + a 3 a.\x\ 


— aia 2 xf) 


Mx,e) 



X2 = T~ f v , (25) 



X 3 

i. 

where 

A(x, e) = det A(x, e) = 1 — e 2 (a 2 a 3 x\ + a 3 aix\ + aia 2 x\) — 2e 3 'a>ia 2 a 3 xix 2 x 3 . (26) 

We will use the abbreviation dET for this map. As always the case for a HK discretization, 
dET is birational, with the reversibility property expressed as / _1 (x, e) = f(x, — e). 

Proposition 12. [E1EIE] The quantities 

_ 1 - e 2 a 3 a l x 2 , _ 1 - e 2 a l a 2 xj _ 1 - e 2 a 2 a 3 x\ 

1 — e A ci\a 2 x\ 1 — e A a 2 a 3 x\ 1 — e A a 3 a\X 2 

are conserved quantities of dET. Of course, there are only two independent integrals since 
F\F 2 F 3 = 1. 

The relation between Fi and the integrals Hi of the continuous time Euler top is 
straightforward: F^ — 1 + e 2 ajifj + 0(e 4 ). As a corollary of Proposition [T21 we find 
that, for any conserved quantity H of the Euler top which is a linear combination of the 
integrals H 1 , H 2 , H 3 , the three functions H/{1 — e 2 ajakX 2 ) are conserved quantities of 
dET. Hereafter (i,j,k) are cyclic permutations of (1,2,3). In particular, the functions 

m = * A ~ akX { (27) 

are conserved quantities of dET. Again, only two of them are independent, since 
aiifi(e) + a 2 H 2 (e) + a 3 H 3 (e) + e A a 1 a 2 a 3 H 1 (e)H 2 (e)H 3 (e) = 0. 
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Proposition 13. [18] The map dET possesses an invariant volume form: 

dx <f)(x) dx\ A d%2 A dx 3 

det 7T = 77~T J u = w ' u = TF~\ > 

ox (p{x) <p{x) 

where <p(x) is any of the functions 

4>{x) = (1 — e 2 aiajX 2 k )(l — e 2 ajakX 2 ) or (1 — e 2 aiO!jX 2 k ) 2 . 

(The ratio of any two functions (f>(x) is an integral of motion, due to Proposition [7J| ). 

The proof is based on formula ([3]) with the matrix A(x, e) given in fT24"j) . Its determinant 
is given in fl26|) . 

A proper discretization of the Wronskian differential equations fl22l) is given by the 
following statement. 

Proposition 14. The following relations hold true for dET: 

' x 2 x 3 - x 2 x 3 = ei/i(e)(xi + Xi), 
x 3 xi - x 3 xi = eH 2 (t)(x 3 + x 3 ), (28) 
xix 2 - xix 2 = eH 3 (e){x 3 + x 3 ) 



1 ■ 



with the functions Hi(e) from (21). 
The proof is based on relations 



2(1 - e 2 a j a k x 2 )(x i + eatXjXk) 

Xi+Xi = J —— r J - , 29 

Afx, e) 



2({(\ r rf n /,./'])(./', • m, ./>/'/,) ^ 
A(x, e) 

which follow easily from the explicit formulas fl2"o]) . They should be compared with 

~_ _ ,~ ~ s _ 2ea i (x j + eoLjXkXi)(x k + eo^XjXj) 

As pointed out in [T7] , a probable way to the discovery of the conserved quantities of 
dET in [8] was through finding the HK bases for this map. In this respect, one has the 
following results. 

Proposition 15. [T7] 

(a) The set $ = (xf, x\, x\, 1) is a HK basis for dET with dimi£$(x) = 2. Therefore, 
any orbit of dET lies on the intersection of two quadrics in M 3 . 

(b) The set $ = (xf, x 2 , xf) is a HK basis for dET with dimif$ (x) = 1. At each 
point x G 1R 3 we have: 

K^ (x) = [ci : c 2 : c 3 ] = [a 2 xl — a 3 x\ : a 3 x\ — 0:1X3 : a\x\ — a 2 x\\. 

Setting c 3 = —1, the following functions are integrals of motion of dET: 

a 3 x\ - a 2 x\ a x x\ - a 3 x\ 
ci(x) = \ f, c 2 (x) = 1 2- ( 32 ) 

OlXg — <X 2 X\ OxXg — ct 2 x{ 

(c) The set $12 = (x 2 , x 2 , 1) is a further HK basis for dET with dimi^$ 12 (x) = 1. At 
each point x G M 3 we have: Kq> l2 (x) = [d± : d 2 : —1], where 

a 2 (l - e 2 a 3 a 1 x\) ai(l - t 2 a 2 a 3 x\) 

di[x) = 5 5 — , d 2 (x) = 5 5 — . (33) 

01X3 — 02^1 a.\X 2 — a 2 xf 
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These functions are integrals of motion of dET independent on the integrals HM) . We 
have: K$>(x) = Kq, Q © K$ 12 . 

Proof. To prove statement (b), we solve the system 

c\x\ + c 2 x\ = x\, 

The solution is given, according to the Cramer's rule, by ratios of determinants of the 
type 



x 2 X 2 



4e(<x,x 2 - a i x'j)(xi + ea 1 x 2 x 3 )(x 2 + ea 2 x 3 x 1 )(x s + ea 3 x 1 x 2 ) 

AHxTe) ' (34) 



(Here we used (15UI) . (1ST]) ). In the ratios of such determinants everything cancels out, 
except for the factors ctjX 2 — a^x 2 , so we end up with The cancelation of the 

denominators A 2 (x,e) is, of course, no wonder, but the cancelation of all the non-even 
factors in the numerators is rather remarkable and miraculous and is not granted by any 
well-understood mechanism. Since the components of the solution do not depend on e, 
we conclude that functions fl32l) are integrals of motion of dET. 
To prove statement (c), we solve the system 

d\x\ + d 2 x\ = 1, 
d\x\ + d 2 x\ = 1. 

The solution is given by eq. fl33|) . due to eq. fl34l) and the similar formula 

,.2 



4eaj(l — e 2 a.ja>kxf)(xi + eaix 2 x 3 )(x 2 + ea 2 x 3 Xi)(x 3 + ea 3 xix 2 ) 

A 2 (x,e) 



1 xf 
1 x\ 

which, in turn, follows from fl29|) and f )30|) . This time the solution does depend on e, but 
consists of manifestly even functions of e. Everything non-even luckily cancels, again. 
Therefore, functions fl33|) are integrals of motion of dET. □ 



Although each one of the HK bases $o ; ^i delivers apparently two integrals of motion 
each pair turns out to be functionally dependent, as 

a.\C\{x) + a 2 c 2 [x) = a 3 , a\di{x) + a 2 d 2 (x) = e 2 aia 2 a 3 . 

However, functions c\, c 2 are independent on d\,d 2 , since the former depend on x 3 , while 
the latter do not. 

Of course, permutational symmetry yields that each of the sets of monomials $23 = 
(x 2 ,, x 2 , 1) and $13 = (xf, x|, 1) is a HK basis, as well, with dim Kq, 23 (x) = dim i^$ 13 (re) = 
1. But we do not obtain additional linearly independent null-spaces, as any two of the 
four found one-dimensional null-spaces span the full null-space K<$>(x). 

Summarizing, we have found a HK basis with a two-dimensional null-space, as well as 
two functionally independent conserved quantities for the HK discretization of the Euler 
top. Both results yield integrability of this discretization, in the sense that its orbits 
are confined to closed curves in IR 3 . Moreover, each such curve is an intersection of two 
quadrics, which in the general position case is an elliptic curve. 
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Proposition 16. Each component xi of any solution of dET satisfies a relation of the 
type Pi(xi,Xi) = 0, where Pi is a biquadratic polynomial whose coefficients are integrals of 
motion of dET: 



Pi(Xi, Xi) = pf ^X-xf + pf\xl + X 2 ) + p^XiXi + pf } = 0, 



with 



pf) = -Ae 2 a 3 a kl pf ] = (l + e^H^e)) (l - e 2 a k H k (e)), 
P? ] = -2(1 - e 2 o^(e)) (l + e 2 a k H k (e)), pf = Ae 2 H 3 (e)H k (e) 
Proof. From eqs. and there follows: 

(Xi Xi 



— ^- + e 2 H?(e)( Xl + Xl f = 2(x)x\ + x)x\). 

It remains to express x 2 and x 2 k through x 2 and integrals Hj(e), H k (e) given in eq. fl27|l . □ 

It follows from Proposition [16] that solutions Xi(t) as functions of the discrete time 
t G 2eZ are given by elliptic functions of order 2 (the order of an elliptic function is the 
number of the zeroes or poles it possesses in a period parallelogram). 

We would like to point out that Propositions [T^andlTolcan be interpreted as existence of 
further HK bases. For instance, according to Proposition!!!] each pair (xjX k —XjX k , Xi+Xi) 
is a HK basis with a 1-dimensional null-space. Similarly, Proposition [16] says that for 
each i — 1, 2, 3, the set x^xl (0 < p, q < 2) is a HK basis with a 1-dimensional null-space. 
Of course, due to the dependence on the shifted variables x, these HK bases consist of 
complicated functions of x rather than of monomials. A further instance of HK bases of 
this sort is given in the following statement. Compared with Proposition [15] it says that 
for dET, for each HK basis consisting of monomials quadratic in x, the corresponding set 
of monomials bilinear in x,x is a HK basis, as well. This seems to be a quite general 
phenomenon, further issues of which will appear later several times. 

Proposition 17. 

(a) The set \I> = (xixi, X2X2, X3X3, 1) is a HK basis for dET with dimi^^(x) = 2. 

(b) The set \l/o = X2X2, is a HK basis for dET with dimi^^ (x) = 1. At 
each point x G M 3 , the homogeneous coordinates q of the null-space K^ (x) = [c\ : C2 : C3] 
are given by 

C{ = (oijxl — a k x 2 ){l — e 2 (oLiOLjX 2 k + a k a{X 2 — aja k x 2 )). 
The quotients Ci/cj are integrals of motion of dET. 

(c) The set \I/ 12 = (xiXi, X2X2, 1) is a further HK basis for dET with dim K^ 12 (x) = 1. 
At each point x G M. 3 , there holds: K^ 12 {x) = [di : d 2 : —1], where 

0^2(1 — e 2 a^a\X2) 1 — e 2 (a.20>zx\ — a^aix 2 + aio^l) 



dx(x) 
d 2 (x) 



ct\x\ — ct2x\ 1 — e 2 {a.20t^x\ + ct 3 ct\x\ — aia 2 x^) ' 
«i(l — e 2 a 2 «3^i) 1 — e 2 (a 3 aixl — a 2 ct3X 2 + aia 2 xl) 



ql\x\ — a.2x\ 1 — e 2 {a.zOL\x\ + c^o^xf — 011012X 2 ) ' 
are integrals of dET. We have: K^(x) = K^ (x) © K^, 12 (x). 
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7. Zhukovski-Volterra SYSTEM 

The gyroscopic Zhukovski-Volterra (ZV) system is a generalization of the Euler top. 
It describes the free motion of a rigid body carrying an asymmetric rotor (gyrostat) [28] . 
Equations of motion of the ZV system read 

xi = aix 2 x 3 + (5 3 x 2 - /3 2 x 3 , 

x 2 = a 2 x 3 xi + (5ix 3 - /3 3 xi, (35) 
x 3 = a 3 x±x 2 + (5 2 xi - /3ix 2 , 

with aii, fa being real parameters of the system. For , /3 2 , /3 3 ) = (0,0,0), the flow (1551) 
reduces to the Euler top (12~T1) . The ZV system is (Liouville and algebraically) integrable 
under the condition 

ai + a 2 + a 3 = 0. (36) 

It can be explicitly integrated in terms of elliptic functions, see [28] and also [3] for a more 
recent exposition. The following quantities are integrals of motion of the ZV system: 

Hi = a 2 x\ - a 3 x\ - 2(p x x\ + (5 2 x 2 + P3X3), 

H 2 = a 3 xj - otixl - 2(fi 1 xi + (5 2 x 2 + $3X3), (37) 

H 3 = olix\ - a 2 x\ - 2(PiXi + (5 2 x 2 + P3X3). 

Clearly, only two of them are functionally independent because of aiHi+a 2 H 2 +a3H 3 = 0. 
Note that 

H 2 - #1 = a 3 C, H 3 -H 2 = aiC, H x -H z = a 2 C, 

with C = x\ + x\ + x\. 

As in the Euler case, the Wronskians of the coordinates Xj admit quadratic expressions 
with coefficients dependent on the integrals of motion: 

x 2 x 3 - x 2 x 3 = H1X1 + xi(/3iXi + f3 2 x 2 + /3 3 x 3 ) + ftC, 
x 3 xi - X3X1 = H 2 x 2 + x 2 (/3iXi + f3 2 x 2 + /3 3 x 3 ) + /3 2 C, (38) 
iix 2 - xii 2 = H 3 x 3 + x 3 ((3ixi + (3 2 x 2 + (3 3 x 3 ) + (3 3 C. 

The HK discretization of the ZV system is: 

xi-xi = eai(x 2 x 3 + x 2 x 3 ) + e(3 3 (x 2 + x 2 ) - e(3 2 (x 3 + x 3 ), 

x 2 - x 2 = ea 2 (x 3 xi + x 3 Xi) + e/3i(x 3 + x 3 ) - e/3 3 (xi +x x ), (39) 

x 3 - x 3 = ea 3 (xix 2 + x{x 2 ) + e(5 2 (xi + x x ) - e/3i(x 2 + x 2 ). 

The map / : x i— > x obtained by solving (1391) for x is given by: 

x = f(x, e) = A _1 (x, e)(l + eB)x, 

with 

(1 -eaix 3 -eaix 2 \ / f3 3 -/3 2 \ 

-ea 2 x 3 1 -ea 2 xi - eB, B = -f3 3 0i \ ■ 
-ea 3 x 2 -ea 3 xi 1 J \ f3 2 -/3 X J 

We will call this map dZV. Formula ([3]) holds true for dZV, as for any HK discretization. 
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7.1. ZV system with two vanishing /^.'s. In the case when two out of three /3fe's 
vanish, say 2 = 3 = 0, the condition fl36|) is not necessary for integrability of the ZV 
system. The functions H 2 and H 3 as given in fl37j) (with f3 2 = (3 3 = 0) are in this case 
conserved quantities without any condition on ct^'s, while their linear combinations H\ 
and C are given by 

1 . 2 2 a 2 + Oi 3 

#i = {a 2 H 2 + a 3 H 3 ) = a 2 x 3 - a 3 x 2 + 2(3% xi, 

Oil Oi 

C=-(H 3 -H 2 )=xl + xl-^±^xl 

Oi Oi 

Wronskian relations ( |38|) are replaced by 

fl «2 + «3 2 , 

£2X3 - x 2 x 3 = H1X1 - pi x\ + PiC, 

Oi 

x 3 x\ - x 3 x\ = H 2 x 2 + (3\Xix 2) 
x\x 2 - x\x 2 = H 3 x 3 + /3iXix 3 . 



(40) 



The HK discretization of the ZV system with (3 2 = (3 3 = turns out to possess two 
conserved quantities (without imposing condition f )36|) ) and an invariant measure. 

Proposition 18. The functions 

a 3 x\ — a>ix\ — 2flxX\ + e 1 0\a\x\ 



1 — e z a 3 ctiX2 

2 2 

H 3 6 



oi\x\ — a 2 x\ — 2/3iXi — e 2 f3la>ix\ 



1 9 2 ' 
1 — e A a\Ct 2 x% > 

are conserved quantities of dZV with (3 2 = (3 3 = 0. 

Proposition 19. The map dZV with (3 2 = j3 3 = possesses an invariant volume form: 

dx (fi(x) dxi A dx 2 A dx 3 

det 7T = 77~T ^ J u = w ' u = JT\ ' 

ox (p[x) <p{x) 

with 4>(x) = (1 — e 2 «3Q;ia;2)(l — e 2 axa 2 x^). 

The conserved quantities of Proposition [IS] appear on the right-hand sides of the fol- 
lowing relations which are the discrete versions of the Wronskian relations (|40j) : 

Proposition 20. The following relations hold true for dZV with (3 2 = (3 3 = 0: 

x 2 x 3 - x 2 x 3 = eci(x\ + xi) + 2ec 2 XiXi + 2ec 3 , 
X3X1 - X3X1 = eH 2 (e)(x 2 + x 2 ) + e0i(xix 2 + x{x 2 ), 
x x x 2 - x x x 2 = eH 3 (e)(x 3 + x 3 ) + e0i(xix 3 + x{x 3 ), 



with 



a 2 H 2 (e) + a 3 H 3 {e) 0t(a 2 + a 3 ) 0i (H 3 (e) - H 2 (e)) 



ci = t , c 2 = - , c 3 



A = 1 + e\a 2 H 3 (e) - f3f) {a 3 H 2 (e) + 2 ). 



Next, we describe the HK bases found in this case. 
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Proposition 21. 

(a) The set $ = (x\, x 2 , x\, x%, 1) is a HK basis for dZV with 02 = 03 = ; with 
dimK<s,(x) = 2. Any orbit of dZV with 02 = 03 = is thus confined to the intersection of 
two quadrics in R 3 . 

(b) The set <3>o = (^1, £ 2 , x 3 , 1) is a HK basis for dZV with (3 2 = 03 = 0, with 
dim K<s> (x) — 1. At each point x G M 3 we have: K^ (x) — [— 1 : d 2 : d 3 : dj], where 

d 2 = (1 - e 2 /3 2 - e 2 a 3 H 2 (e)) , d 3 = (l - e 2 /3 2 + e 2 a 2 H 3 (e)) , 

4 = 7 (H 2 (e)-H 3 (e)). 

(c) The set $23 — (^1) x \-> x ii -0 ^ s a fraszs /or (iZU 02 = 03 = 0, wift 
dim if <j, 23 (x) = 1. At each point 1 £ I 3 we /lave: K$ 23 (x) = [c x : c 2 : c 3 : c 4 ] ; where 

ci = ai(a 3 + e 2 /3 2 a 2 + e 2 a 2 a 3 if 2 (e)) , c 2 = -«i(« 2 + e 2 /3 2 a 3 - e 2 a 2 a 3 if 3 (e)), 
c 3 = -2/3i(a 2 + a 3 ), c 4 = -(a 2 # 2 (e) + a 3 # 3 (e)). 

Unlike the case of dET, we see that here a HK basis with a one dimensional null-space 
already provides more than one independent integral of motion. 
"Bilinear" versions of the above HK bases also exist: 

Proposition 22. The set \I> = (x\Xi, x 2 x 2 , £3X3, x\ + x\, 1) is a HK basis for dZV with 

02 = 03 = 0, with dimi^(x) = 2. The sets 

^0 = (zi^i, x 2 x 2 , X3X3, 1) and ^ 23 = (a; 2 x 2 , x 3 x 3 , x x + Xi, 1) 
are HK bases with one- dimensional null-spaces. 

The following statement is a starting point towards an explicit integration of the map 
dZV with 02 = 03 = in terms of elliptic functions. 

Proposition 23. The component X\ of the solution of the difference equations sat- 
isfies a relation of the type 

P(xi, Xi) = Poxjxl + PiXiXi(x! + Xx) + P2{x\ + x\) + p 3 XiX\ + Pi{x\ + Xl) + £>5 = 0, 

coefficients of the biquadratic polynomial P being conserved quantities of dZV with 02 = 

03 = 0. 

Proof is parallel to that of Proposition [16J 

7.2. ZV system with one vanishing 0^. In the case ^3 = (say) and generic values of 
other parameters, the ZV system has only one integral H 3 and is therefore non-integrable. 
One of the Wronskian relations holds true in this general situation: 

XiX 2 - Xi± 2 = H 3 x 3 + 0\x x x 3 + 2 x 2 x 3 . (41) 

Under condition ( |36|) . the ZV system becomes integrable, with all the results formulated 
in the general case. 

Similarly, the map dZV with 3 = and generic values of other parameters possesses 
one conserved quantity: 

ct\x\ - a 2 x\ - 2{0 x xi + 2 x 2 ) - e 2 (0fai + 0\a 2 )x\ 

n 3 [e) = 2 • 

1 — e z aia 2 x 3 
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Clearly, this fact can be re-formulated as the existence of a HK basis $ = (x\, x 2 ,, x\, x x , x 2 , 1) 
with dimfCj, = 1. The Wronskian relation f T4"Tj) possesses a decent discretization: 

x x x 2 - x x x 2 = eH 3 (e)(x 3 + x 3 ) + e/3x(xix 3 + Xix 3 ) + e(3 2 (x 2 x 3 + x 2 x 3 ). 

However, it seems that the map dZV with j3 3 = does not acquire an additional integral 
of motion under condition (136]) . It might be conjectured that in order to assure the 
integrability of the dZV map with j3 3 = 0, its other parameters have to satisfy some 
relation which is an 0(e)-deformation of (136]) . 

7.3. ZV system with all /3fc's non-vanishing. Numerical experiments indicate non- 
integrability for the map fl39|) with non- vanishing /3k s, even under condition f )36|) . Nev- 
ertheless, some other relation between the parameters might yield integrability. In this 
connection we notice that the map dZV with (oti, a 2 , a 3 ) = (a, —a, 0) admits a polynomial 
conserved quantity 



8. VOLTERRA CHAIN 

8.1. Periodic Volterra chain with N = 3 particles. Equations of motion of the 
periodic Volterra chain with three particles (VC3, for short): 



This system is Liouville and algebraically integrable, with the following two independent 
integrals of motion: 

Hi = x\ + x 2 + x 3l H 2 = xix 2 x 3 . 
There hold the following Wronskian relations: 



Eliminating Xj,Xk from equation of motion for Xi with the help of integrals of motion, 
one arrives at 

x\ = x\[xi - Hi) 2 - AH 2 Xi, 

which yields a solution in terms of elliptic functions. 

The HK discretization of system (T4*2|) (with the time step 2e) is: 

( Xi - Xi = exi(x 2 - x 3 ) + exi(x 2 - x 3 ), 



H = -ax\ - 2{f3 x xi + /3 2 x 2 + /3 3 x 3 ) + e 2 a(/3 2 xi - (5ix 2 ) 2 . 




(42) 



= HiXiXj — 3H 2 . 



(43) 




(44) 



The map / : x 1— > x obtained by solving (JHJ) for x is given by: 

x = f(x, e) = A~ x (x, e)x, 



with 




1 + e(xi - x 3 ) 
ex 3 
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Explicitly: 

~ = x , 1 + 2e( Xj - x k ) + e 2 (( Xj + x k ) 2 - aj) 

1 1 — e 2 {x\ + x\ + 2§ — 2xiX 2 — 2x 2 x 3 — 2x 3 Xi) 
This map will be called CIVC3. From Proposition [2] there follows immediately: 

Proposition 24. The map dVC 3 possesses an invariant volume form: 

dx </>(x) dxi A dx2 A dx^ 
det 7) = 77~T ^ / w = w > w = I7~\ ' 

(72 0(2j 

<p(x) = x\x 2 x 3 . 

Concerning integrability of dVC3, we note first of all that Hi is an obvious conserved 
quantity. The second one is most easily obtained from the following discretization of 
Wronskian relations 



Proposition 25. For the map dVC^, the following relations hold: 

XiXj — XiXj = eHi(xiXj + XiXj) — 6eH 2 (e) (l — ^e 2 Hf) , (46) 

where ^(e) is a conserved quantity, given by 

/ \ X1X2X3 

1 — e 2 {x\ + x\ + x\ — 2xix 2 — 2x 2 x 3 — 2x 3 Xi) ' 

Proof. Define H 2 {e) by eq. (j46]). It is easily computed with explicit formulas fT4"5|) . The 
result given by fT4T|) is a manifestly even function of e and therefore an integral of motion. 

□ 

Proposition 26. 

(a) The set $jj = (xiXj(xi + Xj), x 2 + x 2 , XiXj, Xi + Xj, 1) is a HK basis for the map 
dVC^ with dim K^,^ (x) = 1. In other words, the pairs (xi,Xj) lie on a cubic curve 

P(Xi, Xj) = p X i X j (x i + Xj) + Pi(x 2 + X 2 ) + p 2 XiXj + p 3 (Xi + Xj) + p 4 = 0, 

whose coefficients p m are constant (expressed through integrals of motion). 

(b) The set = (xfx 2 , XjXi{xi + Xi), x 2 + xj, XiXi, x-i + Xi, 1) is a HK basis for the 
map dVC-j with dimi^^-(x) = 1. In other words, the pairs (xi,Xi) lie on a symmetric 
biquadratic curve with constant coefficients (which can be expressed through integrals of 
motion). 

Proof. Statement (a) follows by eliminating x k from fT4"T|) via x k = H\—X{—Xj. Statement 
(b) is obtained with the help of MAPLE; it implies that Xi as functions of t are elliptic 
functions of degree 2 (i.e., with two poles within one parallelogram of periods). □ 

8.2. Periodic Volterra chain with iV = 4 particles. Equations of motion of VC4 are: 

r ±1 = xi(x 2 - x A ), 
x 2 = x 2 (x 3 - Xi), 
x 3 = x 3 (x 4 - x 2 ), 
Xi = 24(21 — 23). 

This system possesses three obvious integrals of motion: H\ = 21+22+23+24, H 2 = xix 3 , 
and H 3 = 2224. One easily finds that x\, x 3 satisfy the differential equation 

x\ = (xj - H lXi + H 2 ) 2 - 4# 3 2?, 



22 MATTEO PETRERA, ANDREAS PFADLER, AND YURI B. SURIS 

while X2, x 4 satisfy a similar equation with if 2 H 3 . This immeadiately leads to solution 
in terms of elliptic functions. There are two types of Wronskian relations: 

X1X3 - X1X3 = 2H 2 (x 2 - X4), X2X4 - X2X4 = 2H 3 (x 3 - £1), (48) 

and 

X1X2 — X\X2 = H\X\X 2 — IH2X2 — 2H 3 X\. 

HK discretization (denoted by dVC 4 ): 

' — _i_ exi(x 2 — x 4 ), 
ex 2 (x 3 -Xx), 
ex 3 (x4 - x 2 ), 



It possesses an obvious integral Hi = x\ + X2 + x 3 + x±. 

Proposition 27. For the map dVC^, the following natural discretization of eqs. ( fjM ) 
holds: 



Xi 


— X\ 


= txi (x 2 


— £4 


x 2 


-x 2 


= ex 2 (x 3 


— X\ 


x 3 


-x 3 


= ex 3 (x4 


- x 2 


X,{ 


— £4 


= ex/Sx\ 


- x 3 



£1X3 - xix 3 = 2eH 2 (e)(x 2 + x 2 - x 4 - x 4 ), 
x 2 x 4 - X2X4 = 2eH 3 (e)(x 1 + x x - x 3 - x 3 ), 



with the conserved quantities 



H,(e) = z * lX3 _ H 3 (e) 



1 — e 2 (x 2 — X4) 2 ' " 1 — e 2 (xi — x 3 ) 2 

Proof. This can be shown directly; the fact that H 2 (e), H 3 (e) are even functions of e, 
assures that they are conserved quantities. One can also show this immediately from 
equations of motion: for instance, multiplying the equations 

Xi Xi x 3 x 3 



1 + e(x2 — X4) 1 — e(x2 — X4) ' 1 — e(x2 — X4) 1 + e(x2 — X4) ' 
shows that H 2 (e) is a conserved quantity. □ 
Proposition 28. 

(a) For the iterates of map dVC4, the pairs (xi,x 2 ) lie on a quartic curve whose coef- 
ficients are constant (expressed through integrals of motion). 

(b) The pairs (xi,Xi) lie on a biquartic curve of genus 1 with constant coefficients 
(which can be expressed through integrals of motion). 

Proof. Statement (a) follows by eliminating x 3 ,X4 from integrals Hi, H 2 (e), H 3 (e). State- 
ment (b) is obtained with the help of MAPLE; it implies that Xi as functions of t are 
elliptic functions of degree 4 (i.e., with four poles within one parallelogram of periods). □ 

Note that the reduction £4 = of the periodic Volterra chain with N = 4 leads to the 
open-end Volterra chain with N = 3 particles. 
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9. Dressing chain (N = 3) 

The three-dimensional dressing chain (DC3, for short) is described by the following 
system of quadratic ordinary equations [27] : 

= x\ - x\ + a 3 - a 2 , 

- x\ — xl + ai — a 3 , (49) 
= x\ — x\ + a 2 — «i, 

with real parameters a^. The system fT4*9|) is (Liouville and algebraically) integrable. The 
following quantities are integrals of motion: 

h = xi + x 2 + £3, 

^2 = + x 2 )(x 2 + x 3 )(x 3 + x{) - a x xi - a 2 x 2 - a 3 x 3 . 

Sometimes it is more convenient to use the following integral instead of I 2 : 

H 2 = x\ + x\ + x\ + 3«iXi + 3a 2 x 2 + 3a 3 x 3 = if — 3I 2 . 

There hold the following Wronskian relations: 

XiXj — Xiij = l\x\ + 2(aj — a>k)Xi + 2(aj — ctk)xj + 3akh — H 2 . (50) 

Excluding Xj,Xk from equations of motion for Xi with the help of integrals of motion, 
one arrives at 

x\ = x\ + Qa 2 x1 + 4a 3 Xj + a 4 , (51) 

with 

a 2 = + Qij + a k - 2ai), a 3 = (aij + a k - a^h + I 2 , 

a 4 = l\ - 2(<x,- + a k )I\ + (atj - a k f - AIJ 2 . 

All three elliptic curves corresponding to ( 15T|) with i = 1,2,3, have equal Weierstrass 
invariants (expressed through the parameters «j and the integrals of motion): 

g 2 = CL4 + 3a 2 , g 3 = a 2 a A — a 2 — a 3 , so that a 3 = — 4a 2 + g 2 a 2 — g 3 . 

The coefficients in (15T|) can be thus parametrized in terms of the Weierstrass elliptic 
function with the invariants g 2 ,g 3 as follows: a 2 = — p(Aj), a 3 = p'(Ai), so that a 4 = 
g 2 — 3p 2 (Ai). One can show that A\ + A 2 + A 3 = (modulo the period lattice), so that 
one can introduce Bi, defined up to a common additive shift, through Ai = Bi — Bi + \. 
The solution of the dressing chain DC3 is then given as 

Xi (t) = ((t - B i+1 ) - C(t - Bi) - C{B t - B l+1 ). 

The HK discretization of system (|49p is: 

%i — x\ = e(x 3 x 3 — x 2 x 2 + a 3 — a 2 ), 

x 2 — x 2 = e(xiXi — x 3 x 3 + a.i — Oi 3 ), (52) 
x 3 - x 3 = e(x 2 x 2 - X1X1 + a 2 - ai). 
The map / : x 1— > x obtained by solving (|52i) for x is given by: 

x = f(x, e) = A _1 (x, e)(x + ec), 
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with 

1 ex 2 — ex 3 
A(x,e)—\-ex 1 1 ex 3 I , c = (a 3 — a 2 , ax — a 3 , a 2 - 



ex\ —ex 2 1 



Explicitly: 



~ _ x i + e ( x t -x] + ack- ocj) + e 2 (hXjX k + (a k - cn^Xj + (acj - aci)x k ) 

1 + e 2 (xix 2 + x 2 x 3 + x 3 x x ) 
This map will be called CIDC3. From Proposition [1] there follows immediately: 
Proposition 29. The map dDC 3 possesses an invariant volume form: 

dx (fi(x) dxi A dx 2 A dx 3 

det IT = TTT ^ / w = w > w = ' 

ax 0(x) <p[x) 

with 4>(x) = 1 + e 2 (x\x 2 + x 2 x 3 + x 3 x\). 

Concerning integrability of dDC3, we note first of all that I± is an obvious conserved 
quantity. The second one is most easily obtained from the following discretization of 
Wronskian relations ( 150]) . 

Proposition 30. For the map dDC 3 , the following relations hold: 



hx k x k + {ai - a k )(xi + Xi) + (aj - a k ){xj + Xj) + 3a k h - H 2 (e), (54) 
where H 2 (e) is a conserved quantity, given by 

H 2 (e) = - 2f H * + * G * (55) 

1 + e^(xix 2 + x 2 x 3 + x 3 xi) 

where 

G 2 = l\x\x 2 x 3 + {pt\X\ + a 2 x 2 + a 3 x 3 ) (x\x 2 + x 2 x 3 + x 3 x{) 

+2Ii(a±x 2 x 3 + a 2 x 3 x\ + 0:3X1X1) 

— {a 2 - a 3 ) 2 x\ - (o 3 - cn) 2 x 2 - («i - a 2 ) 2 x 3 . 

Proof. Define H 2 (e) by eq. (154]) . It is easily computed with explicit formulas (l53j) . The 
result given by (155]) is a manifestly even function of e and therefore an integral of motion. 

□ 

Proposition 31. 

(a) The set $^ = (xf, x|, XiXj(xi + Xj), x 2 , XiXj, x 2 , Xj, Xj, 1) is a HK basis for the 
map dDC 3 with dim K$.. (x) = 1. In other words, the pairs Xi,Xj satisfy equations of 
degree 3, 

Pij{xi, xj) = p Q x\ + pix] + p 2 XiXj(xi +Xj)+ p 3 x 2 + PiXiXj + p 5 x 2 + p 6 Xi + p 7 Xj +p 8 = 0, 

whose coefficients p m = Pm ore constant ( expressed through parameters a k and integrals 
of motion). 

(b) The set = (x™x")mn=o is a HK basis for the map dDC 3 with dim (x) = 1. 
In other words, the pairs Xj,Xj lie on bicubic curves of genus 1: 

3 

Qi{Xi, Xj) ^ ^ Qmn%i 0, 
m,n=0 
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whose coefficients q mn = qmn are constant ( expressed through parameters atk and integrals 
of motion). Moreover, q\ 3 = q 3 ±. 

Proof. Statement (a) follows by eliminating Xk from fl55|) via Xk — I\—X{ — xj. Statement 
(b) is obtained with the help of MAPLE. One can also show that these bicubic curves are 
of genus 1, so that Xi as functions of t are elliptic functions of degree 3 (i.e., with three 
poles within one parallelogram of periods). □ 

10. Coupled Euler tops 

In [6j a remarkable mechanical system was introduced, which can be interpreted as a 
chain of coupled three-dimensional Euler tops. The differential equations governing the 
system are given by: 

x\ = oix 2 x 3 , 

X2j = <X3j-l%2j-l%2j+l, 

%2j+l = Cl2,jX2jX2j-l + a 3 j + iX2j+2 x 2j+3, (56) 
^2^+2 — C^3j+2 x 2j+l x 2j+3j 
&2N+1 — OtZN%2NX2N-l, 

with real parameters Oj. Each triple of variables (a^j-i, %2j, %2j+i) can be considered as 
a 3D Euler top, coupled with the neighboring triple (X2.7+1, x 2j+2, ^+3) via the variable 
%2j+i- We will denote system (156]) by CETjy. It has N + 1 independent conserved 
quantities: 

Hi = ct2x\ — a.\x\, 

Hj = a^^a^^xlj^ - a 3j _4a 3j _ia;2 i _ 1 + a 3j _ 4 a 3j ^ 2 xlj, 2 < j < N, 

Hn+1 = &3NX2N ~ °i3N-l x 2N+l- 

Nothing is known about the possible Hamiltonian formulation of this system, and therefore 
about its integrability in the Liouville-Arnold sense. 

For N = 1 system (1561) reduces to the usual Euler top (12T1) . We will consider in detail 
the HK discretization of the system CET2 given by 

( ±1 = a 1 x 2 x 3 , 
x 2 = a 2 x 3 xi, 

x 3 = a 3 xiX2 + 0:4X4X5, (57) 
x 4 = o 5 x 5 x 3 , 
X5 = 05X3X4. 

It can be interpreted as two Euler tops, described by the two sets of variables (xi,X2,X3) 
and (x3,X4,X5), respectively, coupled via the variable X3. It has three independent inte- 
grals of motion: 

Hi = a 2 x\ — OJ1X2, H 3 = o 6 xf — 05X5, 
H 2 = 0305X2 — o 2 o 5 X3 + o 2 o 4 X4, 

and it can be solved in terms of elliptic functions. We will be mainly interested in its 
particular case which is superintegrable. 
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Proposition 32. If the coefficients on satisfy the following condition, 

Oi0 2 = o 5 o 6 , 

then the system CET 2 is superintegrable: it has two additional integrals, 

i?4 = 05X2X5 — O2X1X4, H 5 = 0:5X1X5 — 01X2X4, 

and among the functions Hi, . . . , H5 there are four independent ones. 
In this case, the variable 23 satisfies the following differential equation: 



.1. o 



+ — ) ( 

V at 2 at 5 J V 



Q'l 



030:4 



01O2X3 H -H 2 + a 3 Hi - a A H 3 - 

o 5 J a 2 a 5 

so that its time evolution is described by an elliptic function of degree 2. 
The HK discretization of the system CET2 reads: 

( - " eoi(x 2 x 3 + x 2 x 3 ), 

eo 2 (x 3 xi + X3X1), 

eo 3 (xiX2 + x\x 2 ) + £04(^4X5 + X4X5), 
eo 5 (x 5 x 3 + X5X3), 
eo 6 (x 3 x 4 + X3X4). 
The map / : x 1— > x obtained by solving (1601) for x is given by: 

A' 1 ' 



H 



Xi 


— Xi 


x 2 


- x 2 


X3 


- x 3 


X4 


— X4 


x 5 


- x 5 



X 



x,e x, 



with 



A{x,e) 



( 1 eoiX3 eoix 2 

£02X3 1 602X1 

eo 3 x 2 eo 3 Xi 1 

605X5 

V 



eo 6 x 4 eo 6 x 3 
This map will be called dCET 2 in the sequel. 







£04X5 €04X4 

1 eo 5 x 3 
1 



Proposition 33. The functions 
HAe 



ol 2 x\ 



a\x\ 



H»(e) 



a 6 xl 



a 5 x 2 5 



e 2 a 5 a e xl ' 



(5f 



(59) 



(60) 



1 — e 2 0i0 2 X3 " 1 

are conserved quantities of the map dCET 2 . 

Proposition [M] gives only two independent integrals of motion for the map dCET 2 . 
Numerical experiments indicate that the third integral does not exist in general. The 
situation is different under condition ( |58|) . Note that in this case the denominators of the 
integrals Hi(e) and H 3 (e) coincide. 

Proposition 34. If condition l[58\) holds, then the map dCET 2 has in addition to Hi(e) 
and H 3 (e) also the following conserved quantities 



H 4 (e) 



# 2 (e) 
05X2X5 • 



03O5XI — 02O5X3 + 02O4XI 



e 2 0i02x| 



02X1X4 



e 2 0i02x| 



H,(e) 



O5X1X5 — O1X2X4 
1 — e 2 o 1 o 2 x 2 ; 
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There are four independent functions among i?i(e), . . . , H^{e). 
We now present HK bases for the map 0ICET2. 



Proposition 35. Under condition §58$ , the map dCET 2 has the following HK bases. 

(a) The set $ = (x 2 , x\, x§, x 2 , x\, 1) is a HK basis with dimi£$(x) = 3. 

(b) The sets $1 = (xf, x 2 ,, x\, 1) and $2 = (^3, ^4, xf, 1) are HK bases with one- 
dimensional null-spaces. At each point x G IR 5 we have: K^^x) = [e\ : e 2 '■ e 2 a\a 2 : —1] 
and K$ 2 (x) = [e 2 0i0 2 '■ fi '■ fs '■ — 1]- The functions and f\ are conserved quantities 
given by 

o 2 (l — e 2 a>ia 2 x\) Oi(l — e 2 a\a 2 x^) 

C\ o 9 , 62 o ^ , 

a 2 xf — a>ix 2 ct 2 x\ — ct\X 2 

and 

_ o 5 (l - e 2 aia 2 x\) a^l - e 2 a 1 a 2 xl) 

J 4 9 2 ' i5 o 2 ' 

0:5X4 — 04X5 5 X| — O4X5 

The set i/if fraszs wrat/i a one-dimensional null-space. At each 

point x & M. 5 we have: Kq> 3 (x) = [g\ : g 2 : o 5 : — o 4 ]. T/ie functions gi are conserved 
quantities given by 

a 3 a 5 x^ — a 2 a 5 xl + a 2 a^x\ O3O5X 2 — 0105X3 + O1O4X4 

9i 2 2 ' 92 o 2 ■ 

O2X4 — O1X2 2 X^ — OiXg 

Similar claim hold for the sets (x 2 , x 2 , x|, xf), (xf, x|, xf, x 2 ), and (x 2 , x|, x|, xf). 

(c) T/ie sei \1/ = (xi, x 2 , x 3 , x 4 , x 5 ) is a HK basis with dimK^(x) = 2. 

(d) The sets = (xi, X2, X4) ; ^ 2 = {x\, x 2 , X5) are HK bases with one- dimensional 
null-spaces. At each point x G M 5 we have: K^^x) — [c\ : c 2 — 1], K^ 2 (x) = [di : d 2 : 
— 1]. The functions ci,c 2 and di,d 2 are conserved quantities given by 

02X1X4 — 05X2X5 05X1X5 — 01X2X4 

C\ K 9 ' o 2 ' 

O2XJ — O1X2 O2XJ — OiXg 

while 

, _ ct 2 _ Oi 

c(i — — c 2 , a 2 — — c%. 

o 5 o 5 

A similar claim holds for the sets (x2,X4,X5) and {x\^x^x^). 

We see that the map dCET 2 under condition (JSSJ) possesses four functionally indepen- 
dent conserved quantities. It might look paradoxical that 3> U \IMs a HK basis with a 
5-dimensional null-space, thus imposing seemingly 5 restrictions on any orbit of the map, 
which would yield O-dimensional invariant sets (instead of invariant curves in the continu- 
ous time case). The resolution of this paradox is that the five restrictions are functionally 
dependent on their common set (i.e., along any orbit). In other words, the HK basis $U\I/ 
is not regular. This is the first and the only instance of a non-regular HK basis in this 
paper. 

The map dCET 2 possesses, in its (super)-integrable regime, an invariant volume form: 



Proposition 36. Under condition §oo\) , the map dCET 2 preserves the following volume 
form: 

dx 0(x) dx\ A dx 2 A dx% A g?X4 A dx$ 

det 7T = TT\ j w = w, u = -7— , 

ox <p{x) <p{x) 
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with <p(x) = (1 — e 2 aia 2 x 3 1 ) 3 . 

In this regime, the solutions can be found in terms of elliptic functions, as the following 
statement shows. 



Proposition 37. Under condition foty) . the component x 3 of any orbit of the map dCET 2 
satisfies a relation of the type 

Q(x 3 , x 3 ) = q xlxl + gix 3 x 3 (x 3 + x 3 ) + q 2 (xl + x\) + q 3 x 3 x 3 + + x 3 ) + q 5 = 0, 

coefficients of the biquadratic polynomial Q being conserved quantities of dCET 2 . 

This statement is a proper discretization of eq. (|59p . 

11. Three wave system 
The three wave interaction system of ordinary differential equations is [2]: 

Z\ = \OL\Z 2 Z 3 , 

z 2 = \a 2 z 3 z u (61) 

Z 3 = \OL 3 Z X Z 2 . 

Here z = (zx, z 2 , z 3 ) G C 3 , while the parameters a, of the system are supposed to be real 
numbers. If (i, j, k) stands for any cyclic permutation of (123), then we can write system 
(I6"T|) in the abbreviated form 

^ = iaiZjZ k , (62) 
Writing Zi = Xi + iy iy i = 1,2, 3, we put system (1621 into the form 

x% oci^Xjyh -\- yjX^, 

i \ (63) 

Ui = ai(XjX k - yjy k )- 

System (163]) is completely integrable and can be integrated in terms of elliptic functions. 
It has three independent integrals of motion: quadratic ones, 

Hi = aj\z k \ 2 - a k \zj\ 2 , 

among which there are only two independent ones because of ot\Hi + a 2 H 2 + a 3 H 3 = 0, 
and a cubic one, 

K = - (z x z 2 z 3 + zxz 2 z 3 ) = y\z{z x z 2 z 3 ). 



The HK discretization of system (l62|) reads 

Zi — Zi = ieai(zjZ k + ZjZ k ), 

or, in the real variables (a^,^), 

Xi-Xi = eai(xjy k + x d y k + yjX k + yjX k ), 
Vi-y% = eai(xjX k + xjx k - y^y k - yjy k ). 

In the matrix form, this can be put as 

A(x,y,e)(z\ = Pj & ft) = f(x,y,e) = A~ l (x } y } e) P 
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where 



/ 


1 


-eaiy 3 


-tct x y 2 





—£0:1X3 


—eaix 2 \ 




-ea 2 y 3 


1 


-ea 2 yi 


—ea 2 x 3 





—ea 2 Xi 




-ea 3 y 2 


-ea 3 yi 


1 


—ea 3 x 2 
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The birational map / : MP — > M 6 will be called d3W hereafter. 

Proposition 38. The map d3W has three independent conserved quantities, namely, any 
two of 

Hi(e) 

supplied with any one of 

9te(ziz 2 Z3)(l ~ e 2 afcai|zj| 2 )(l - e 2 a^|z fc | 2 ) 

K A e ) = T7 — — -v , 

A(z,z,e) 

where 

A(z,z,e) = det A(x, y, e) = 1 — 2e 2 (a 2 a 3 \zi\ 2 + a 3 ai\z 2 \ 2 + aia 2 \z 3 \ 2 ) 

+ e (a^o^l^il +«3Q;i|-2 2 | + aia^l^l ) — Ae a 1 a 2 a 3 \zi\ \z 2 \ \z 3 \ . 
Proposition 39. The map d3W possesses an invariant volume form: 

dz (f)(z) dxi A dx 2 A dx 3 A dy l A dy 2 A dy 3 

det d- z = W) * /w=w ' u= W) ' 

where <f>(z) = A(z, z, e) . 

Next, we give the results on the HK bases for the map d3W, which yield a complete 
set of integrals of motion. 

Proposition 40. 

(a) The sets $j = (|zj| 2 , \zk\ 2 , 1), % = 1,2,3, are HK bases for the map d3W with 
dim K^.(z) = 1. At each point z e C 3 there holds: K^.(z) = \d\ : d 2 : —1], where the 
coefficients 

di(z) = "^""l ' , d 2 (z) 
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are integrals of motion of the map d3W. They are functionally dependent because of 
ajdi(z) + a k d 2 (z) = e 2 a\a 2 a 3 . 

(b) The sets = (D\e(z 1 z 2 z 3 ) , \zi\ 2 , 1), % = 1,2,3, are HK bases for the map d3W with 
dim Ky^z) = 1. At each point z G C 3 there holds: Ky^z) — [e± : e 2 : — 1], where the 
coefficients 

A(z,z,e) 



e±(z) = 

^fKz(ziZ 2 z 3 ){l - e 2 (-ajCtklzil 2 + a k ai\zA 2 + aiaAz k \ 2 ) ) 

e 2 (s) = 



(l - e 2 (-a^afcl^l 2 + afca^] 2 + a t a j\z k \ 2 ) ) 
are independent integrals of motion of the map d3W. 
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12. Lagrange top 

Lagrange top was the second integrable system, after Euler top, to which the HK 
discretization was successfully applied [12J. We reproduce and re-derive here the results 
of that paper, and add some new results. 

Equations of motion of the Lagrange top are of the general Kirchhoff type: 

m = mx V m H + px V P H, 
p = px V m H, 

where m = (mi, rri2, m 3 ) T and p = (pi,P2, P3) T - Any Kirchhoff type system is Hamiltonian 
with the Hamilton function H = H(m,p) with respect to the Lie-Poisson bracket on c(3)*, 

{mi^rrij} = e ijk m k , {m u Pj} = UjkPk, {Pi,Pj} = 0, 

and admits the Hamilton function H and the Casimir functions 

Ci = Pi + Pi + pI, C 2 = m x px + m 2 p 2 + m 3 p 3 , (65) 

as integrals of motion. For the complete integrability of a Kirchhoff type system, it should 
admit a fourth independent integral of motion. 

The Hamilton function of the Lagrange top (LT) is H = H\/2, where 

Hi = ml + ml + ami + 2^p 3 . 

Thus, equations of motion of LT read 

' rhi = (a - l)m 2 m 3 + 7p 2 , 
rh 2 = (1 — «)mim 3 - Tpi, 

™ 3 = °' (66) 
p\ = ap 2 m 3 -p 3 m 2 , 

P2 = P3^i - aptm 3 , 
p 3 = p x m 2 -p 2 m 1 . 

It follows immediately that the fourth integral of motion is simply H 2 = m 3 . Traditionally, 
the explicit integration of the LT in terms of elliptic functions starts with the following 
observation: the component p 3 of the solution satisfies the differential equation 

Pi = PM, (67) 

with a cubic polynomial P 3 whose coefficients are expressed through integrals of motion: 

P 3 (p 3 ) = (Hi - ami - 2 7 p 3 )(Ci - p 2 3 ) - (C 2 - m 3 p 3 f . 

We mention also the following Wronskian relation which follows easily from equations 
of motion: 

(riiipi - mipi) + (rn 2 p 2 - m 2 p 2 ) + (2a - l)(rn 3 p 3 - m 3 p 3 ) = 0. (68) 
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Applying the HK discretization scheme to eqs. fl66|) . we obtain the following discrete 
system: 

mi - mi = e(a - l)(m 2 m 3 + m 2 m 3 ) + e^(p 2 + P2), 
m 2 - m 2 = e(l - a)(mim 3 + mim 3 ) - ey{pi + pi), 
fh 3 -m 3 = 0, 

pi - Pi = ea(p 2 m 3 + p 2 m 3 ) - e(p 3 fh 2 + p 3 m 2 ), 
p 2 -p 2 = e(p 3 mi + p 3 mi) - ea(pim 3 + pim 3 ), 
P3 - P3 = e(Pi^2 + Pim 2 - p 2 mi - p 2 mi). 
As usual, this can be solved for (m,p), thus yielding the reversible and birational map 



141 = f(x, e) 
I 



A 1 (x, e)(l + eB)x, where x = (mi, m 2 , m 3 ,pi,p 2 ,p 3 ) 



A(x,e) 
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This map will be called dLT in the sequel. Obviously, m 3 serves as a conserved quantity 
for dLT. The remaining three conserved quantities can be found with the help of the HK 
bases approach. A simple conserved quantity can be found from the following statement 
which serves as a natural discretization of the Wronskian relation 



Proposition 41. The set V = {rhipi — mipi, m 2 p 2 — m 2 p 2 , rh 3 p 3 — m 3 p 3 ) is a HK basis 
for the map dLT with dimi^r(^) = 1- At each point x € R 6 we have: Kr(x) = [1:1: 63], 
where b 3 is a conserved quantity of dLT given by 



(2a — l)m 3 + e 2 (a — l)m 3 {m\ + m 2 ) + e 2 7(mipi + m 2 p 2 ) 



where 



Ai 



m 3 Ai 



1 + e 2 a(l — a)m\ — e 2 jp 3 . 



(69) 
(70) 



Proof. A straightforward computation with MAPLE of the quantity 

(m^i - mipi) + (m 2 p 2 - m 2 p 2 ) 
(m 3 p 3 - m 3 p 3 ) 

leads to the value (1691) . It is an even function of e and therefore a conserved quantity. □ 



Further integrals of motion were found by Hirota and Kimura. We reproduce here their 
results with new simplified proofs. 
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Proposition 42. [12] 

(a) The set $ = (m\ + m 2 , p\m\ + ^2^2, Pi +p\-, p\, Ps, 1) is a HK basis for the map 
dLT with dimK^(x) = 3. 

(b) The set $i = (1, p 3 , p 2 , m\ + m 2 .) is a HK basis for the map dLT with a one- 
dimensional null-space. At each point x G M 6 we have: K^ >1 (x) = [co : C\ : C2 : — 1]. The 
functions Cq,Ci,C2 are conserved quantities of the map dLT, given by 

m 2 + m 2 + 2 7 p 3 + e 2 4 4) + e 4 c { 6) + e 6 ^ + e 8 c^ 0) 
Co — 



c-i 



c 2 



AiA 2 

2 7 (l - e 2 a(l - a)m 2 3 ) (l + e 2 cf + e 4 c { 2 4) + e 6 cf ] ) 



AiA 2 

;V(l + e 2 4 2) +6 4 4 4) + e 6 4 6) ) 
AiA 2 



Here A 1 is given in (7§, and A 2 = 1 + e 2 A { 2 2) + e 4 A^ 4) + e 6 A^ 6) ; coefficients A^ and c[ g) 
are polynomials of degree q in the phase variables. In particular: 

cf 1 = ml + m 2 2 + (1 - 2a + 2a 2 )m 2 - 27p 3 , 
Ag 2 '* = ml + ml + (1 — 3a + 3a 2 )m 2 — 7p 3 . 

(c) The set $2 = (1, P3, P3, m\P\ + ^2^2) a HK basis for the map dLT with a one- 
dimensional null-space. At each point x G IR 6 we have: Kq, 2 (x) = [do '■ d\ : 0*2 : — 1]- The 
functions 0*0, di, 0*2 are conserved quantities of the map dLT, given by 

_ mipi + m 2 p2 + m 3 p 3 + e 2 df ) + e 4 dg 6) + e 6 dg 8) + e 8 4 10) 
do — 



d 1 
d 2 



AiA 2 

m 3 + e 2 af } + e 4 af ] + e 6 dj 7) + e 8 af } 
A X A 2 

e 2 7 (l - a)m 3 (l + e 2 cf + e 4 4 4) + e 6 c 2 6) ) 
A X A 2 



where are polynomials of degree q in the phase variables. In particular, 

= 7(mipi + ni2P2) — 7 (3 — 2a)m 3 p 3 + am 3 (m 2 + m 2 ) + (1 — 3a + 3a 2 )m 3 . 

(d) The set $ 3 = (1, p 3 , p 2 , p\ + p 2 ,) a HK basis for the map dLT with a one- 
dimensional null-space. At each point x G M 6 we have: K<f> 3 (x) — [eo : e\ : e% : — 1]. T/ie 
functions eo,e\,e2 are conserved quantities of the map dLT, given by 

Pi + Vl + Vl + e 2 e^ 4) + e 4 e< 6) + e 6 e 8) + e 8 e { W) 
eo — 



AiA 2 

2e 2 (ei 3) + 6 2 ef ) +eM 7) +e 6 ef ) ) 
A X A 2 

(1 + e 2 (l - a) 2 m|) (l + e 2 ^ 2) + e 4 4 4) + e 6 cf) 



AiA 
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where ejj. are polynomials of degree q in the phase variables. In particular, 

ef = lip\ +pl+ pf) ~ (1 - a)m 3 (mipi + m 2 p 2 + m 3 p 3 ). 

Proof, (b) We consider a linear system of equations 

(c + cxp 3 + c 2 pl) of(m,p,e) = (m 2 + m 2 ) o f l (m,p, e), (71) 

for all ieZ. Numerically one sees that it admits a unique solution, and one can identify 
the linear relation 

\ie 2 ci= (l-e 2 a(l-a)ml)c 2 . (72) 

The system of three equations for three unknowns cq, ci, C2 consisting of flTT]) with z = 0, 1 
and f!72|) can easily be solved with MAPLE. Its solutions are even functions of e, which 
proves that they are integrals of motion. 

(c) This time we consider the linear system of equations 

(d + dip 3 + d 2 pl) o f(m,p, e) = (m x pi + m 2 p 2 ) o f(m,p, e), (73) 

for all i £ Z. Numerically we see that it admits a unique solution, and we can identify 
the linear relation 

jd 2 = (1 - a)m 3 c 2 . (74) 



The system of three equations for the three unknowns do,di,d 2 consisting of (J73J) for 
i — 0, 1 and of fT74l) with C2 already found in part b) can easily be solved with MAPLE. 
Its solutions are even functions of e and therefore are integrals. 

(d) Completely analogous to the last two proofs: we solve the linear system of three 
equations for the three unknowns e , e±, e 2 , consisting of the equations 

(e + eip 3 + e 2 p 2 3 ) o f(m,p,e) = (pj+pj) o f(m,p,e), 

for i = 0,1, and of the linear relation 

e 2 7 2 e 2 = (l + e 2 (l-a) 2 m 2 )c 2 , 

and verify that they are even functions of e. □ 

We note that for a = 1 the integrals do, d\, d 2 simplify to 

mipx + m 2 p 2 + m 3 p 3 m 3 + e 2 ^(m 1 p l + m 2 p 2 ) 

d = = , di = , d 2 = 0. 

1 - e 2 ^p 3 1 - e^7p 3 

It is possible to find a further simple, in fact polynomial, integral for the map dLT. 
Proposition 43. [12] The function 

F = m 2 + m 2 , + 2^p 3 - e 2 ((l - a)m 3 mi + 7Pi) 2 - e 2 ((l - a)m 3 m 2 + 7p 2 ) 2 , 
is a conserved quantity for the map dLT. 

Proof. Setting C = 1 — e 2 (l — a) 2 m\, D = —2e 2/ ~f(l — a)m 3 , E = — e 2 7 2 , one can check 
that Cci + Ddi + Ee\ = and Cc 2 + Dd 2 + Ee 2 = —2j. This yields for the conserved 
quantity F = Ccq + Ddo + Eeo the expression given in the Proposition. □ 
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Considering the leading terms of the power expansions in e, one sees immediately that 
the integrals cq, do, eo, and m 3 are functionally independent. Using exact evaluation 
of gradients we can also verify independence of other sets of integrals. It turns out 
that for a ^ 1 each one of the quadruples {do, d\, d 2 , m^) and {eo, e±, e 2 , ^3} consists of 
independent integrals. 

A direct "bilinearization" of the HK bases of Proposition H2] provides us with an alter- 
native source of integrals of motion: 

Proposition 44. The set 

\J = ( mi mi + m 2 m 2 , Pimi + V\ m \ + P2m 2 + p 2 m 2 , P1P1 + P2P2, P3P3, P3 + P~3, 1) 
is a HK basis for the map dLT with dimi^^(x) = 3. Each of the following subsets of ^/ , 
^1 = (1, p 3 + p 3 , p 3 p 3 , mifht + m 2 rh 2 ), 
^2 = (1, P3 + P3, P3P3, miPi + rhiPi + m 2 p 2 + m 2 p 2 ), 
^3 = (1, P3 + P3, P3P3, P1P1 + P2P2), 
is a HK basis with a one-dimensional null-space. 

Concerning solutions of dLT as functions of the (discrete) time t, the crucial result is 
given in the following statement which should be considered as the proper discretization 
ofeq. flEU). 

Proposition 45. [12] The component p% of the solution of difference equations sat- 
isfies a relation of the type 

Q(P3,Ps) = QopIpI + QiPsMps + Ps) + <?2<>3 + Ps) + QsPsPs + Qa(P3 + Ps) + Q5 = 0, 

coefficients of the biquadratic polynomial Q being conserved quantities of dLT. Hence, 
Pzif) is an elliptic function of degree 2. 

Although it remains unknown whether the map dLT admits an invariant Poisson struc- 
ture, we have the following statement. 

Proposition 46. The map dLT possesses an invariant volume form: 

dx (p(x) dm 1 A dm 2 A dm 3 A dpi A dp 2 A dp 3 

det 7) = TT\ ^ / u = u i w = -7— , 

ox <p{x) <p{x) 

with 4>(x) = A 2 (x,e). 



13. KlRCHHOFF CASE OF THE RIGID BODY MOTION IN AN IDEAL FLUID 

The motion of a rigid body in an ideal fluid is described by Kirchhoff equations f )M|) 
with H being a quadratic form in m = (mi,m 2 ,m 3 ) T G M 3 and p = (pi,p 2 ,p 3 ) T G M 3 . 
The physical meaning of m is the total angular momentum, whereas p represents the total 
linear momentum of the system. A detailed introduction to the general context of rigid 
body dynamics and its mathematical foundations can be found in |14j . 

The integrable case of this system found in the original paper by Kirchhoff [13] and 
carrying his name is characterized by the Hamilton function H = Hi/ 2, where 



Ht = ai(m? + ml) + a 3 m\ + b x (p\ + p\) + b 3 p 2 3 . 
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The differential equations of the Kirchhoff case are: 

rni = (a 3 - a 1 )m 2 m 3 + (6 3 - h)p 2 p 3 , 



m 2 = (ai - a 3 )mim 3 + (61 - b 3 )p x p 3 , 
m 3 = 0, 

Pi = a 3 p 2 m 3 - a 1 p 3 m 2 , 
P2 = a 1 p 3 m 1 - a 3 pim 3 , 
P3 = a 1 (p 1 m 2 -p 2 m 1 ). 



(75) 



Along with the Hamilton function H and the Casimir functions ( 165]) . it possesses the 
obvious fourth integral, due to the rotational symmetry of the system: H 2 = m 3 . Tradi- 
tionally, the explicit integration of the Kirchhoff case in terms of elliptic functions starts 
with the following observation: the component p 3 of the solution satisfies the differential 
equation 

Pi = Pa(Pb), 

with a quartic polynomial P4 whose coefficients are expressed through integrals of motion: 

Plipz) = ai (Hi - a 3 m 2 3 - b^d - p§) - b 3 p\) (d - p\) - a\(C 2 - m 3 p 3 ) 2 . 

We mention also the following Wronskian relation which follows easily from equations 
of motion: 

ai(ihipi - mxpi) + ai(m 2 p 2 - m 2 p 2 ) + (2a 3 - ai)(rn 3 p 3 - m 3 p 3 ) = 0. (76) 

Applying the HK approach to f!75|) . we obtain the following system of equations: 

mi - mi = e(a 3 - ai)(m 2 m 3 + m 2 m 3 ) + e(6 3 - h) (p 2 p 3 + p 2 p 3 ) , 
m 2 -m 2 = e(ai - a 3 )(m 1 m 3 + mj^) + e(&i - & 3 )(piP 3 + PiP 3 ), 
m 3 -m 3 = 0, 

Pi - Pi = ta 3 {p 2 m 3 + ^2^3) - eai{p 3 m 2 + p 3 m 2 ), 
p 2 -p 2 = eai(p 3 mi +p 3 mi) - ea 3 (pim 3 + pim, 3 ), 
P3 -P3 = ea 1 (p 1 m 2 +p 1 m 2 ) - eai(p 2 m x +p 2 m 1 ). 

As usual, these equations define a birational map x = f(x, e), x = (m,p) T . We will refer 
to this map as dK. Like in the case of dLT, m 3 is a conserved quantity of dK. A further 
"simple" conserved quantity can be found from the following natural discretization of the 
Wronskian relation f[7"6l). 



Proposition 47. The set T = (rhipi — mxpi, m 2 p 2 — m 2 p 2 , rh 3 p 3 — m 3 p 3 ) is a HK basis 
for the map dK with dimKr(x) = 1. At each point x G M 6 we have: Ky(x) = [1:1: —73], 
where 7 3 is a conserved quantity of dK given by 

73 = (77) 
aiAi 



where 



A = ai - 2a 3 + e 2 ai( a i - a 3 )(m\ + mj) + e 2 a 1 a 3 (b 1 - b 3 )(p\ + pj), (78) 
Ai = 1 + e 2 a 3 (ai - a 3 )m\ + e 2 a 1 (b 1 - b 3 )p\. (79) 
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Proof. Like in the case of dLT, we let MAPLE compute the quantity 

(fhipi - mipi) + (m 2 p 2 - rn 2 p 2 ) 

73 = -pz , 

{m 3 p 3 - m 3 p 3 ) 

which results in (1771) . an even function of e and therefore a conserved quantity. □ 
Interestingly enough, this same integral may also be obtained from another HK basis: 

Proposition 48. The set $o — (wf + m\,p\ + p\,p\, 1) is a HK Basis for the map dK 

with dim K<s> (x) = 1. The linear combination of these functions vanishing along the orbits 
can be put as A — 7301 Ai = 0. 

Proof. The statement of the Proposition deals with the solution of a linear system of 
equations consisting of 

(c 1 (ml + m 2 2 ) + c 2 {p 2 1 +p 2 2 ) + c 3 pl)of i (m,p,e) = 1 (80) 

for all i 6 Z. We solve this system with i = —1,0, 1 (numerically or symbolically), and 
observe that the solutions satisfy a 3 {bi — b 3 )c\ = ai(ai — a 3 )c 2 . Then, we consider the 
system of three equations for c±, c 2 , c 3 consisting of the latter linear relation between c±, 
c 2 , and of eqs. (I8"U|) for i — 0, 1. This system is easily solved symbolically (by MAPLE), 
its unique solution can be put as in the Proposition. Its components are manifestly even 
functions of e, thus conserved quantities. □ 

Proposition 49. 

(a) The set $ = (m\ + m 2 ,, p\m x + p 2 m 2 , p\ + p 2 , pi, p 3 , 1) is a HK basis for the map 
dK with dim.K^(x) = 3. 

(b) The set $1 = (1, p 3 , p\, m\ + m 2 ) is a HK basis for the map dK with a one- 
dimensional null-space. At each point x G M 6 we have: K^^x) = [co : C\ : c 2 : —1]. The 
functions co,ci,c 2 are conserved quantities of the map dK, given by 

ai {m\ + ml) - (h - b 3 )pl + e 2 4 4) + e^cf + e 6 4 8) + e 8 c 10) 
Co — 



Cl 



c 2 



aiA!A 2 

2e 2 a 3 (b 1 - b 3 )m 3 (C 2 + e 2 cf + e 4 c[ 6) + e^cf) 



AiA 2 

(h - 63) (1 + e 2 cf + e 4 c 2 4) + e 6 cj> 6) + e 8 cf) 
aiAiA 2 



where Ai is given in ( f7P| ), and A 2 = 1 + e 2 A^ +e 4 A 2 A ' } +e 6 Af \ coefficients cj^ and A^ 
are homogeneous polynomials of degree q in the phase variables. In particular: 

c 2 2) = -2a\(m\ + ml) - (a\ - 2ai<2 3 + 2>al)m\ + ai(&i - b 3 )(p\ + pi) - a^bi - b 3 )pl, 
A^ = a\(m\ + ml) + (a\ - 3aia 3 + 3al)ml - ai(bi - b 3 ) (p\ + pi) + a x (bi - 6 2 )pg. 

(c) The set $ 2 = (l,p 3 ,pl, m\p\ + m 2 p 2 ) is a HK basis for the map dK with a one- 
dimensional null-space. At each point x G M 6 we have: K<$> 2 (x) = [do : d\ : d 2 : —1]. The 
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functions do,di,d,2 are conserved quantities of the map dK, given by 

? a + e 2 4 4) + e 4 af + e 6 af + e*d£ 0) 



d 
d\ 
d2 



AiA 2 

m 3 ( - 1 + gag + e 4 af } + e 6 af > + e 8 af } ) 
AiA 2 

01(63 - Qi)e 2 (C7 2 + e 2 cj 4) + eM 6) + e 6 cf } ) 
AiA 2 



where d k are homogeneous polynomials of degree q in the phase variables. In particular, 
df ] = -aia 3 (m 2 + m 2 .) - (a 2 - 3aia 3 + 3a 2 )m 2 + {a x - 03) (61 - 63) (Pi +P2) _ 3a i ( & i ~ ^3)^3- 

(d) Tae sei $3 = (1, P3, pf, p 2 + p 2 ) a 6aszs /or i/ie map dK with a one- 
dimensional null-space. At each point x G M 6 we have: K$ z {x) = [eo : e\ : e 2 : — 1]. 
The functions eo,ei,ez are conserved quantities of the map dK, given by 

d + e 2 ei 4) + + + e^ 0) 
e o — 



ei 



A X A 2 

2e 2 a!(a 3 - a 1 )m 3 (C 2 + e 2 cf + e 4 ^ + e 6 cj 8) ) 



AiA 

_l +e2e ( 2 ) +e 4 e (4). 

e 2 



_1 + e 2 e ( 2 ) + e 4 e W + ,6,(6) + e8e (8) 



AiA 2 

where are polynomials of degree q in the phase variables. In particular, 



,( 2 ) 



-a\(m\ + ml) - (2a 2 - Aa x a ?J + 3a 2 )m 2 + 2ai(&i - 63) (Pi +P 2 ) - a i(^i - ^3^3- 



Proof. Statement (b) is proven using direct calculation. Statements (c) and (d) then 
follow analogously to Proposition H2] from the existence of linear relations between c\ and 
d 2 , as well as between c\ and e%. □ 

One can show that each of the sets {co,ci,c 2 }, {d , d\, d 2 }, and {eo,ei,e 2 } consists 
of three independent integrals of motion. Moreover, each of the sets {co, Ci, c 2 , UI3} and 
{eo, ei, e 2 , m 3 } consists of four independent integrals. As further important results, me 
mention that Propositions HH (on the "bilinear" HK bases), H5] (on the invariant bi- 
quadratic curve for (^3,^3)), andH6] (on the invariant measure) hold literally true for the 
map dK. 

14. Clebsch case of the rigid body motion in an ideal fluid 

Another famous integrable case of the Kirchhoff equations was discovered by Clebsch [3] 
and is characterized by the Hamilton function H = Hi/2, where 



Hi = (m, Am) + (p, Bp) = - ^2(a k m 2 k + b k p 2 k ), 

k=l 
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where A 



diag(ai, a 2 , 03) and B 

h- 



= diag(&i,& 2 ,&3 
h h - b 3 63 



0. 

a 3 a\ a 2 

This condition is also equivalent to saying that the quantity 

h 



satisfy the condition 



6 



(81) 



2) 



CLi(aj - a k ) 

takes one and the same value for all permutations (i, j, k) of the indices (1,2,3). 

For an embedding of this system into the modern theory of integrable systems see 
[161 [22] . Note that the Kirchhoff case (ai = a 2 and 61 = b 2 ) can be considered as a 
particular case of the Clebsch case, but is special in many respects (the symmetry resulting 
in the existence of the Noether integral m 3 , solvability in elliptic functions, in contrast to 
the general Clebsch system being solvable in terms of theta-functions of genus 2, etc.). 
Equations of motion of the Clebsch case are: 

m = m x Am + p x Bp , 
p = p x Am, 



3) 



or in components 



rhi = (a 3 - a 2 )m 2 m 3 + (b 3 
m 2 = (ai - a 3 )m 3 mi + (bi 
rii3 = (0-2 - ai)mim 2 + (b 2 
p 1 = a 3 m 3 p 2 - a 2 m 2 p 3 , 
p 2 = axmipz - a 3 m 3 pi, 
p 3 = a 2 m 2 pi - a x m x p 2 . 



b 2 )p 2 p 3 , 
h)P3Pu 



Condition (18 ip can be resolved for (X<i clS 

b 2 - b 3 



a 1 



a 2 



bi 



a 3 



h - b 2 

UJi - UJ 2 



UJ 2 — UJ 3 U 3 — UJ\ 

For fixed values of Ui and varying values of bi, equations of motion of the Clebsch case 
share the integrals of motion: the Casimirs C\, C 2 , cf. eq. flo^]) . and the Hamiltonians 



pU 



m. 



UJi — OJk 



+ 



mt 



Ui — UJi 



There are four independent functions among Ci, Ii, because of C\ = I\ + I 2 + I 3 . Note 
that H\ = b\Ii + b 2 I 2 + 63/3. One can denote all models with the same UJi as a hierarchy, 
single flows of which are characterized by the parameters bi. Usually, one denotes as "the 
first flow" of this hierarchy the one corresponding to the choice bi = UJi, so that Oj = 1. 
Thus, the first flow is characterized by the value 9 = 00 of the constant 



14.1. First flow of the Clebsch system. The first flow of the Clebsch hierarchy is 
generated by the Hamilton function H = Hi/ 2, where 



Hi 



m\ + m\ + m\ + Uipj + UJ 2 p\ + uo 3 p\. 



The corresponding equations of motion read: 

m = p x Qp, 
p = p x m, 
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where Q = diag(o;i, 002, 003) is the matrix of parameters, or in components: 

' mi = (u 3 - w 2 )p 2 p 3 , 

m 2 = (ui - w 3 )p3Pi, 

m 3 = (w 2 - uJx)pip 2 , 

Pi = m 3 p 2 - m 2 p 3 , 

p 2 = mip 3 - m 3 pi, 
k p 3 = m 2 pi - mip 2 . 

The fourth independent quadratic integral can be chosen as 

H 2 = Uxrnl + u 2 ml + uj 3 m\ - u 2 0J 3 p\ — ^3^iPl ~ ^i^2P 3 - 

Note that H\ = uili + w 2 / 2 + cu 3 / 3 , Hi = —oj 2 uj 3 Ii — U3U1I2 — wio; 2 / 3 . 
We mention the following Wronskian relation: 

(rhxpx - mxpx) + (m 2 p 2 - m 2 p 2 ) + (m 3 p 3 - m 3 p 3 ) = 0, (85) 

which holds true for the first Clebsch flow. 

The HK discretization of the first Clebsch flow (proposed in [21]) is: 

' m 1 -m 1 = e(u 3 - w 2 )(p 2 p 3 + p 2 p 3 ), 
m 2 - m 2 = e(ui - u 3 ){p 3 pi + P3P1), 
fh 3 - m 3 = t(u 2 - tO\){p\P2 + P1P2), 
Pi - Pi = e(rh 3 p 2 + m 3 p 2 ) - e(m 2 p 3 + m 2 p 3 ), 
P2~P2 = e(mip 3 + mip 3 ) - e(m 3 pi + m 3 pi), 

, P3~P3 = e(rh 2 pi + m 2 pi) - e(mip 2 + mxp 2 ). 

As usual, it leads to a reversible birational map x = f(x,e), x = (m,p) T , given by 
f(x,e) = A~ 1 (x,e)x with 



/ 


1 











ew 23 p 3 


e^ 23 p 2 \ 







1 





ew 3 ip 3 





eu 3 ip! 










1 


ew i2 p 2 


euj 12 pi 














1 


-em 3 


em 2 









epi 


em 3 


1 


—em! 


\ 


eP2 


-epi 





-em 2 


emi 


1 / 



where the abbreviation UJij — UJ{ — OJj IS used. This map will be referred to as dC. 

A "simple" conserved quantity can be found from the following natural discretization 
of the Wronskian relation (1851) . 

Proposition 50. The set T = {rhipi — mipi, m 2 p 2 — m 2 p 2 , rh 3 p 3 — m 3 p 3 ) is a HK basis 
for the map dC with dimi^r^) = 1. At each point x G M 6 we have: K-p(x) = [ei : e 2 : e 3 ], 
where 

ei = 1 + e 2 (ui - <jjj)pj + e 2 (uii - u k )p 2 k . (86) 

The conserved quantities Cijcj can be put as tijtj = (1 + e 2 WjJ)/(l + e 2 ujjJ), where J is 
a nice and symmetric integral, 

j = pI+pI+ pI 

1 - e 2 (wipf + UJ 2 p\ + u 3 pl) ' 
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Remarkably, it can be obtained also from a different (monomial) HK basis, see part b) of 
the following statement. 

Proposition 51. [UJ 

(a) The set of functions $ = {p\, p\, p\, m\, m\, mj, m\p\, m 2 p 2 , m 3 p 3 , 1) is a HK 
basis for the map dC with dim K^(m,p) = 4. Thus, any orbit of the map dC lies on an 
intersection of four quadrics in M 6 . 

(b) The set of functions $ — (Pi? v\i p\-, 1) a HK basis for the map dC with 
dim Kq> (m,p) = 1. At each point {m,p) e R 6 there holds: 

K& (m,p) = [ei : e 2 : e 3 : -(pI+pI+pI)] 

1 o "^2 2 

— + e ui : — + e'u 2 : -j + e u 3 : — 



with the quantities given in 
(c) The sets of functions 

$1 = (Pi, pIj P3' m i> m i ""4 m iPi), (87) 

$2 = (p?, P 2 , PiL m l> m 2, m 3> m 2P2), (88) 

$3 = (Pi, Pi, pI, m\, ml, ml, m 3 p 3 ), (89) 

are HK bases for the map dC with dim K^ (m, p) = dim K$ 2 (m,p) = dim K$ 3 (m,p) = 1. 
eac/i point (m,p) e M 6 i/iere holds: 

K$ 1 (m,p) = [«i : a 2 : «3 : «4 : «5 : «6 : — 1], 

K^ 3 (m,p) = [71 : 72 : 73 : 74 : 75 : 76 : -1], 

where ctj,/3j, and jj are rational functions of (m,p), even with respect to e. They are 
conserved quantities of the map dCS. For j = 1, 2, 3 ; they are of the form 

_ + e 2 h^ + e 4 h^ + e 6 h^ + e s h^ + e 10 h^ 
~ 2e2(pf+p|+p|)A ' 

where h stands for any of the functions ctj, (3j, jj, j = 1, 2, 3, 

A = mipi + m 2 p 2 + m 3 p 3 + e 2 A (4) + e 4 A (6) + e 6 A (8) , 

and the corresponding h^ 2q > , A^ are homogeneous polynomials in phase variables of 
degree 2q. For instance, 

a x — Oi — 1\, a 2 — a 3 — —ii, 



df ] = ~h, 0? = Ci - h, /f ) = -h, 

(2) r ( 2 ) r ( 2 ) ^ r 

7i = -h, 72 = - J 3, % = Ci - I 3 . 



For j = 4, 5, 6, the functions atj, (3j, jj are given by 



/ a 4 


«5 


a 6 \ 




04 


05 


06 


H 


\74 


75 


76 J 





D Ai/[u)i-u 3 ) A 1 /(ui-u 2 y 

A 2 /(u 2 -u 3 ) D A 2 /(u 2 -u 1 ) 

A 3 /(u 3 -u 2 ) A 3 /(uj 3 -uji) D 
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where 



1 + e'Af + e 4 A£> + e s A^ + e 8 ^' 
2e 2 A 

Pi + p\ + P% + e 2 ^ (4) + e 4 £> (6) + e 6 D 
2A 



| (4) 



(6) 



.8/1(8) 



(8) 



are homogeneous polynomials of degree 2q in phase variables, for instance, 

A^ 2) = m 2 + m 2 , + m 2 + (u 2 + uj 3 - 2u 1 )pl + (lo 3 + oj x ~ 2w 2 )p 2 + ( w i + w 2 - 2w 3 )p 2 . 

The four conserved quantities J, ci\, fi\ and 71 are functionally independent. 

Our paper [T7j contains a much more detailed information about the HK bases of 
the map dC, for instance, the further basis with a one- dimensional null-space: O = 
{pi, p 2 , p 3 , mipi, rri2P2, ^3^3)- However, the following finding about the "bilinear" ver- 
sions of the above bases is new. 



(90) 
(91) 
(92) 
(93) 



Proposition 52. Each of the sets of functions 
*o = (PiPu P2P2, P3P3, 1), 

^1 = (pijPi, P2P2, P3P3, m 1 m 1 , m 2 m 2 , m 3 m 3 , m x p x + miPi), 
^2 = (P1P1, P2P2, P3P3, mi mi, m 2 m 2 , m 3 m 3 , m 2 p 2 + m 2 p 2 ), 
*3 = (piPi, P2P2, P3P3, mpi, m 2 m 2 , m 3 m 3 , m 3 p 3 + m 3 p 3 ), 
zs a i/i^ frasis /or the map dC with a one- dimensional null-space. 

14.2. General flow of the Clebsch system. The HK discretization of the flow fl83l 

reads _ _ 

m — m = e{m x Am + m x Am +px 5p + px 5p ) , 

p — p = e(p x Am +px Am) , 

in components: 

( mi - mi = e(a 3 - a 2 )(m 2 m 3 + m 2 m 3 ) + e(6 3 - b 2 )(p 2 p 3 + P2P3), 
m 2 -m 2 = e(d! - a 3 )(m 3 mi + m 3 mi) + e(6i - & 3 )(P3Pi + P3P1), 
fn 3 -m 3 = e(a 2 - ai)(m l m 2 + mim 2 ) + e(6 2 - &i)(pip 2 +P1P2), 
p\—Pi = ea 3 (m 3 p 2 + m 3 p 2 ) - ea 2 (m 2 p 3 + m 2 p 3 ), 
p 2 -p 2 = eai(mip 3 + mip 3 ) - ea 3 (m 3 pi + m 3 pi), 
^ P3~P3 = ea 2 (m 2 pi + m 2 Pi) - ea x {rh x p 2 + mip 2 )- 



(94) 



In what follows, we will use the abbreviations 6^- = 
system (194)) defines an explicit, birational map / : 



b,; — bj and a,-,- = a,- — a,-. The linear 



f(m,p,e) 



M~ 



lm,p,e) 



where 



M(m,p, e) 



/ 1 


ea 23 m 3 


ea 23 m 2 





tb 23 p 3 


eb 23 p 2 \ 


ea 3 im 3 


1 


ea 3 imi 


eb 31 p 3 





e&3iPi 


ea 12 m 2 


eai 2 mi 


1 


eb 12 p 2 


tb 12 p! 








ea 2 p 3 


-ea 3 p 2 


1 


-ea 3 m 3 


ea 2 m 2 


-eaip 3 





ea 3 pi 


ea 3 m 3 


1 


— eaimi 


V eaiP2 


-ea 2 pi 





-ea 2 m 2 




1 / 
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This map will be denoted dGC in what follows. 

A "simple" integral of the map dGC can be obtained by discretizing the following 
Wronskian relation with constant coefficients, which holds for the general flow of the 
Clebsch system flMl) : 

Ai(mxpi - miji) + A 2 (m 2 p 2 - m 2 p 2 ) + A 3 (m 3 p 3 - m 3 p 3 ) = 0, 

with 

Ai = d^dj -\- diCLfc — CLjCLfc. 

Proposition 53. The setT = (rhipi—mipi, m 2 p 2 — 1TL2P2, ^3P3 — ^3^3) is a HK basis for 
the map dGC with dimi^r(^) = 1- At each point x G 1R 6 we have: Kr(x) = [ei : e2 : 63], 
where, for k) = c.p.(l, 2, 3), 

ei = A { + e 2 ai(bi - bj)A k Qj + e 2 ai(bi - b k )AjO k , 

with 

r\ 2 1 a i 2 

UdjCL k 

(Recall that 9 is defined by eq. (f&g)): we assume here that 9 7^ 00.) 

As in the case of the first flow, the integrals e^/ej can be expressed through one sym- 
metric integral: e%jej = (Ai — 9aiL)/(Aj — 9ajL), where 

_ a 2 a 3 AiQi + a 3 aiA 2 Q2 + aia 2 A 3 Q 3 
l + e 2 9a 1 a 2 a 3 (& l + Q 2 + 3 ) 

The quantities and the integral L can be also obtained from a different (monomial) HK 
basis, given in part b) of the following Proposition. 

Proposition 54. 

(a) The set $ = (pf, p 2 ,, mf, m 2 , m\, m\P\, m2p 2 , m 3 p 3 , 1) is a HK basis for the 
map dGC with dim K<p(m,p) = 4. Thus, any orbit of the map dGC lies on an intersection 
of four quadrics in R 6 . 

(b) The set of functions $ — (Pi, v\i Pi> m ii m 2i m i> 1) i s a ^K basis for the map 
dGC with dim Kq, Q (m, p) — 1. At each point (m,p) G M 6 there holds: 

K$ (m,p) = [a 2 a 3 ei : a 3 a x e 2 : a x a 2 e 3 : (a 1 /9)e 1 : {a 2 /9)e 2 : {a 3 /9)e 3 : -e ], 

if/iere 

e = a 2 a 3 Ai6i + a 3 aiv4 2 6 2 + aia 2 A 3 Q 3 

is an integral of motion of the continuous time flow (84\ ). 

(c) The sets of functions [Ety- fifty) are HK bases for the map dGC with one-dimensional 
null-spaces. 

(d) Each of the sets of functions ^ = (piPi, p 2 p 2 , P3P3, mimi, m 2 m 2 , m 3 m 3 , 1) and 
< fP7)) -< fP^) zs a /Of fraszs /or i/ie map d(7C wrat/i a one- dimensional null-space. 
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15. SU(2) RATIONAL GAUDIN SYSTEM WITH N = 2 SPINS 

The Gaudin system [7] describes an interaction of N quantum spins yi, i = 1, . . . , N, 
with a homogeneous constant external field p. Its classical version is given by the following 
quadratic system of differential equations [19] : 

1 < i < N, (95) 

where G su(2) ~l 3 ,p£ su(2) ~ R 3 is a constant vector, and pairwise distinct numbers 
Xi are parameters of the model. The flow (|95p is Hamiltonian with respect to the Lie- 
Poisson bracket of the direct sum of TV copies of su(2), admits 2N independent conserved 
quantities in involution: the N Casimir functions 

Ck = (yk,Vk), 

and the following N Hamiltonians: 

where (•, •) denotes the scalar product in su(2) ~ M 3 . Note that the Hamilton function of 
the flow ([15]) is 

N ^ N N 

k=l i,j=l i=l 

In [151119] it has been proved that a contraction of N simple poles to one pole of order TV 
provides the integrable flow of the so called one-body rational su(2) tower, whose simplest 
instance with N = 2 describes the dynamics of the three-dimensional Lagrange top in the 
rest frame. 

We consider here the HK discretization of thes flow (l9"5"j) with N = 2. We set y x = 
(xi, x 2 , x 3 ) T , y 2 = (zi, z 2 , z 3 ) T , and choose the constant gravity vector p = (0, 0, 1) T . We 
thus obtain the following system of differential equations: 





= x 2 z 3 


- X 3 Z 2 + XiX 2 


x 2 


= X 3 Z! 


— X\Z 3 — Ai^i 


x 3 


= X\Z 2 


— X 2 Z\, 


Zl 


= z 2 x 3 


— Z 3 X 2 + \ 2 Z 2l 


z 2 


= Z 3 Xx 


— z\x 3 — X 2 zi, 


Z3 


= ZiX 2 


- z 2 x u 



(96) 



with Ai, A2 being real parameters. The system (I96p has the following four independent 
integrals of motion: 

C\ = X2 ~\~ X2 ~\~ x^, C 2 = z^ + Z2 ~\~ z^ , 

X X Zi + X 2 Z 2 + X 3 Z 3 X1Z1 + X 2 Z 2 + X 3 Z 3 

Hi = x 3 H , H 2 = z 3 H . 

Al — A 2 A 2 — Ai 
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0. 



(97) 



Note that the quantity Hi + H 2 = x 3 + z 3 is a linear integral of motion. We mention also 
the following Wronskian relation with constant coefficients: 

(x 3 + Z 3 ){x X Z X - XxZx + x 2 z 2 - x 2 z 2 ) + (Ai + \ 2 +x 3 + z 3 )(x 3 z 3 - x 3 z 3 ) 

The HK discretization of (!96|) reads: 

xi — x\ = e(x 2 z 3 + x 2 z 3 - x 3 z 2 - x 3 z 2 ) + e\i(x 2 + x 2 ), 
x 2 - x 2 = e(x 3 zi + x 3 zi - xiz 3 - x{z 3 ) - eAi(xi + xi), 
x 3 - x 3 = e(xiz 2 + xiz 2 - x 2 z\ - x 2 zi), 
Zx — z\ = e(z 2 x 3 + x 2 x 3 - z 3 x 2 - z 3 x 2 ) + t\ 2 {z 2 + z 2 ), 
z 2 - z 2 = e(z 3 X! + 23X1 - Z!X 3 - zix 3 ) - e\ 2 (Ii + zi), 
z 3 - z 3 = e(zix 2 + z\x 2 - z 2 x\ - z 2 xi). 

The map / : x 1— > x obtained by solving fl98l) for x is given by: 

x = f(x, e) = A^(x, e)(l + eB)x, 
where x = (x\, x 2 , x 3 , Zi, z 2 , z 3 ) T , and 



(9* 



A(x,e) 
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-ezi 
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ex 2 


— ex 1 










ez 3 


-ez 2 
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-ex 3 


ea;2 




-ez 3 





ez 1 




ex 3 
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— exi 
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-ex 2 
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-A 2 
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0/ 





This map will be called dG in the sequel. The quantity x 3 + z 3 is obviously preserved by 
the map dG. Other conserved quantities may now be found using the HK bases approach. 
A "simple" integral follows, as in the previous sections, by discretizing the Wronskian 
relation ([9T 



Proposition 55. The set of functions F = ^X\Z\ —x{zi, x 2 z 2 —x 2 z 2 , x 3 z 3 — x 3 z 3 ) is a HK 
basis for the map dG with the one- dimensional null-space K-p(x, z) = [x 3 + z 3 : x 3 + z 3 : I], 
with 

, , _ Ai + A 2 + x 3 + z 3 + e 2 \ x {x\ + x 2 2 ) + e 2 A 2 (^ + z\) + e 2 (Ai + \ 2 ){x\Z\ + x 2 z 2 ) 
{X ' Z) ~ l-e 2 (A l3 ;3 + A 2 ^3 + AiA 2 ) • 

A full set of integrals is found in the following Proposition. The roles of the variables 
Xi and Zi are not quite symmetric there, and interchanging them is of course admissible 
but does not lead to new integrals of motion. 

Proposition 56. 

(a) The set $ = {x\ + x 2 , z\ + z 2 , z 3 , X\Z\ + x 2 ^ 2 , z 3 , 1) is a HK basis for the map dG 
with dimfCj>(x) = 3. 
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(b) The set $1 = (1, z 3 , z 2 , x\ + x?,) is a HK basis for the map dG with a one- 
dimensional null-space. At each point (x,z) G M 6 we have: K® 1 (x,z) = [cq : C\ : c 2 : — 1]. 
The functions co,ci,c 2 are conserved quantities of the map dG, given by 



c = 



Cl 



c 2 



+ xi - 2x3^3 - 4 + + gM. + gff; + is*. 

AiA 2 

2(x 3 + z 3 ) (l + e 2 cS 3) + e 4 cS 5) + e 6 cf } + e 8 cf } ) 
A^ ' 

(l + 6 2 4 2) + 6 4 4 4) + 6 6 cf + 6 8 4 10) ) 
A^ ' 



where 



A 1 = 1 - e 2 ^ + A 2 z 3 + AxA 2 ), A 2 = 1 + e 2 A< 2) + e 4 A< 4) + e 6 A< 6) . 
i/ere and are polynomials of degree q in the phase variables. In particular: 
c 2 2) = -{x\ + x 2 2 + xf) - (z\ + zl + z\) -2(x 1 zi + x 2 z 2 + x 3 z 3 ) -2(X 2 x 3 + Aiz 3 ) - (A 2 + A 2 ), 

and = c 2 2) + Aix 3 + X 2 z 3 + A X A 2 . 

(c) The set $ 2 = (1, z 3 , z\, X\Z\ + x 2 z 2 ) is a HK basis for the map dG with a one- 
dimensional null-space. At each point (x,z) G M 6 we have: K$ 2 (x,z) = [d : d\ : d 2 : — 1]. 
The functions d ,d 1 ,d 2 are conserved quantities of the map dG, given by 

, x lZl + x 2 z 2 + x 3 z 3 + (A 2 - X 1 )z 3 + e 2 4 4) + ^ + ^ + ^4° } 
do = 



d 1 = 



AiA 2 

Ai - A 2 - x 3 - z 3 + e 2 df + t A df + 6 6 rff } + 6 8 4 9) 

AiA 2 



1 + 6 2 4 2) + 6 4 4 4) + 6 6 cf + 6 8 4 8) 

d2 ~ A^ ' 

where dfj^ are polynomials of degree q in the phase variables. 

(d) The set <E> 3 = (1, z 3 , z 2 , z\ + z 2 ) is a HK basis for the map dG with a one- 
dimensional null-space. At each point (x,z) G IR 6 we have: K$ 3 (x,z) = [eo : e\ : e 2 : —1]. 
The functions eo,e±,e 2 are conserved quantities of the map dG, given by 

4 + 4 + 4 + ^ + e 4 e< 6) + e Q ef + e 8 e { w) 
60 " ' 



2e 2 ( eS 2) + e 2 eS 4) + e*ef + e«ef ] 



ei = 



AiA 

1 + e 2 4 2) + 

e 2 = 



1 + e 2 e ( 2 ) + e 4 e (4) + £ 6 e (6) + e 8 e (8) 



A X A 2 

where e^ are polynomials of degree q in the phase variables. 
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It can be shown that each of the sets {co, ci, C2, 23 + z 3 }, {do, d\, d 2 , £3 + £3} and 
{eo, ei, e 2 , %3 + £3} contains four independent integrals. 

Analogously to the situation for the map dLT, it is possible to obtain a polynomial 
integral and an invariant volume form for the map dG. 

Proposition 57. The function 

11 e 2 
G = -(xi + zi) 2 + -(x 2 + z 2 ) 2 + A1X3 + X 2 z 3 - — ((A1X1 + \ 2 zi) 2 + (Aix 2 + A 2 2 2 ) 2 ), 

is a conserved quantity for the map dG. 

Proposition 58. The map dG possesses an invariant volume form: 

, d(x, z) 6(x, z) .„ dx\ A dx 2 A da; 3 A dz\ A rf^ A dz 3 

det 7- = — f oj = u), u - 



d(x, z) (j)(x, z) <f>{x, z) 

where <p(x, z)) = A 2 (x, z). 

Explicit integration of the map dG could be based on the following claim. 

Proposition 59. The component £3 of the solution of the dG map satisfies a relation of 
the type 

Q{x 3 , x 3 ) = q xlxl + qix 3 x 3 (x 3 + x 3 ) + + xf) + q 3 x 3 x 3 + q 4 (x 3 + x 3 ) + q 5 = 0, 

coefficients of the biquadratic polynomial Q being conserved quantities of dG. Thus, x 3 (t) 
is an elliptic function of degree 2. An analogous statement holds for the component z 3 . 

16. Conclusions 

The initial motivation for this study was the hope that HK discretization would pre- 
serve the integrability for all algebraically integrable systems. This was formulated as 
a conjecture in [17J. The list of integrable discretizations given in the present overview 
contains more than a dozen issues and is rather impressive. It includes systems integrable 
in terms of elliptic functions as well as those integrable in terms of theta-functions of 
genus g = 2 (Clebsch system). This list might look like a convincing argument in favor 
of the integrability conjecture. However, at present we have also found examples which 
indicate that this conjecture might be wrong (e.g., the Zhukovski-Volterra system with 
all /3k 7^ 0, or integrable chains, Volterra and dressing ones, with a big number of parti- 
cles, say N > 5). We do not have rigorous proofs of the non- integrability in these cases, 
but the numerical evidence is rather strong. Since HK discretizations are (probably) not 
always integrable, the big number of integrable cases becomes still more intriguing: it 
is hard to imagine that all their common features come as a pure coincidence. We are 
after a theory which would clarify the problem of integrability of HK discretizations, but 
at present such a theory seems to remain a rather remote goal. Still, even without a 
general framework, HK discretizations represent a new fascinating chapter in the theory 
of integrable systems: we are now in a possession of a big and potentially growing stock 
of birational maps, integrable in terms of Abelian functions, highly non-trivial from the 
point of view of algebraic geometry and very different in nature from anything known be- 
fore. The immediate goal of this review will be achieved if HK discretizations will attract 
attention of experts in the theory of integrable systems and in algebraic geometry. 
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